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ABSTRACT 

We address the issue of when generalized quantum dynamics, which 
is a classical symplectic dynamics for noncommuting operator phase space 
variables based on a graded total trace Hamiltonian H, reduces to Heisen- 
berg picture complex quantum mechanics. We begin by showing that when 
H = Trif , with H a Weyl ordered operator Hamiltonian, then the gen- 
eralized quantum dynamics operator equations of motion agree with those 
obtained from H in the Heisenberg picture by using canonical commuta- 
tion relations. The remainder of the paper is devoted to a study of how 
an effective canonical algebra can arise, without this condition simply being 
imposed by fiat on the operator initial values. We first show that for any 
total trace Hamiltonian which involves no noncommutative constants, there 
is a conserved anti-self-adjoint operator C with a structure which is closely 
related to the canonical commutator algebra. We study the canonical trans- 
formations of generalized quantum dynamics, and show that C is a canonical 
invariant, as is the operator phase space volume element. The latter result 
is a generalization of Liouville's theorem, and permits the application of sta- 
tistical mechanical methods to determine the canonical ensemble governing 
the equilibrium distribution of operator initial values. We give arguments 
based on a Ward identity analogous to the equipartition theorem of clas- 
sical statistical mechanics, suggesting that statistical ensemble averages of 
Weyl ordered polynomials in the operator phase space variables correspond 
to the Wightman functions of a unitary complex quantum mechanics, with 
a conserved operator Hamiltonian and with the standard canonical commu- 
tation relations obeyed by Weyl ordered operator strings. Thus there is a 
well-defined sense in which complex quantum field theory can emerge as a 
statistical approximation to an underlying generalized quantum dynamics. 
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1. Introduction and Brief Review of 
Generalized Quantum Dynamics 

In several recent publications one of us (S.L.A.) formulated [1, 2] and with collabo- 
rators elaborated [3, 4] an operator dynamics, called generalized quantum dynamics, which 
gives a symplectic dynamics for general noncommutative degrees of freedom. This per- 
mits the direct derivation of equations of motion for field operators, dispensing with the 
conventional canonical procedure of "quantizing" a classical theory. Although we observed 
that generalized quantum dynamics in a complex Hilbert space is compatible with canonical 
quantization, the precise connection between the two formalisms was not established, and it 
is this issue which we address in the present paper. We will not in fact find it necessary to 
restrict ourselves to complex Hilbert space, and the derivations and conclusions given here 
apply (with some specific differences which we discuss) in quaternionic Hilbert space and 
real Hilbert space as well. 

Generalized quantum dynamics can be given in either Lagrangian or Hamiltonian 
form, and for brevity we review only the Hamiltonian formalism, since this is what we will 
need. We shall assume an underlying Hilbert space which is the direct sum of bosonic and 
fermionic subspaces, and a grading operator (— 1) F with eigenvalue 1(— 1) for states in the 
bosonic (fermionic) subspace. For a general operator O, we define the graded trace operation 
TrC by 

TrO =ReTr(-l) F = Re^(n|(-l) F C|n) 

(1) 

=ReJ2( n \°\ n ) -Re^(n|C|n) , 

n,B n,F 

with the subscripts B, F on the sums indicating summations over bosonic and fermionic 
states, respectively. We call operators bosonic if they commute with (— 1) F and fermionic if 
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they anticommute with (— 1) F . Given sufficient convergence, it is then easy to see that TrO 
vanishes if O is fermionic, and that Tr obeys the cyclic property 



with the +(— ) sign holding when 0(\) and 0(2) are both bosonic (fermionic). 

Let now {q r (t)}, {p r (t)}, r = 1,...,N be a set of operator phase space variables. 
For each r, we assume that q r and p r are either both bosonic or both fermionic, but we 
make no a priori assumption about the commutativity of the phase space variables with one 
another. Letting A[{g r }, {p r }} be a polynomial (or Laurent expandable) operator function 
of the phase space variables, we define the real number-valued total trace functional A by 



Although noncommutativity of the phase space variables prevents us from simply differenti- 
ating A with respect to them, we can use the cyclic property of Tr to define derivatives of 
A by forming SA and cyclically reordering all the operator variations 5q r ,5p r to the right, 
giving the fundamental definition 



in which 5A/5q r and 5A/5p r are themselves operators. 

Let us now introduce an operator Hamiltonian H[{q r }, {p r }] and a corresponding 
total trace Hamiltonian H = TrH, which generates the dynamics of the phase space variables 
via the operator Hamilton equations 



TrO (1) (2) = ±TrO (2) (1) 



(2) 



A[{q r },{p r }] = TrA[{q r .},{ Pr }} 



(3) 




(4) 



with e r = 1(— 1) according to whether q r and p r are bosonic (fermionic) , and with the dot 

denoting the time derivative. (As discussed in Refs. [1, 2], this Hamiltonian formulation can 

be derived from a total trace Lagrangian action principle, in strict analogy with standard 

derivations of classical mechanics.) If A[{g r }, {p r }, t] is an arbitrary total trace functional 

which can have an explicit time dependence, as well as an implicit time dependence through 

its dependence on the phase space variables, then a simple application of Eqs. (1-5) shows 

that the total time derivative of A is given by 

dA dA r A . , 

1F = lr + {A.H} , (6a) 

where we have denoted by {A, B} the generalized Poisson bracket defined by 

{A B}-TVVe f 5A5B 5B5A\ 

Since the generalized Poisson bracket is antisymmetric in its arguments, by taking A to 
be the total trace Hamiltonian H, which has no explicit time dependence, we learn from 
Eq. (6a) that H is a constant of the motion. 

In addition to its antisymmetry, the generalized Poisson bracket can also be shown 
[3] to satisfy the Jacobi identity 

= {A,{B,C}} + {C,{A,B}} + {B,{C,A}} . (7) 

As a consequence, the phase space flows in generalized quantum dynamics exhibit many 
features [4] analogous to those of ordinary classical mechanics. In order to exhibit the 
symplectic structure of the generalized Poisson bracket, we use the cyclic property of Tr to 
rewrite Eq. (6b) as 



{A, B} = Tr 



^(5A5B_5A5B\ / SA 5B 5A5B\ 

^ \ <5<?r $Pr $Pr Sq r ) \ Sq r 5p r 5p r 5q r J 



(8) 



with the subscripts B, F on the sums respectively indicating summations over the bosonic 
and fermionic phase space variables. If we now introduce the notation x\ — q±, x-i — p±, x% — 
92,^4 = Vii •••,^2A r -i = Qn,X2n = Pn for the operator phase space variables, Eq. (8) can be 
compactly rewritten as 

2iV 

and similarly, the operator Hamilton equations of Eq. (5) can be compactly rewritten as 

2N 

in-., 

' 5x 



x r = ^uj rs — . (9ft) 

s=l 



(Henceforth, we will not explicitly indicate the range of summation indices; the index r on 
q r , p r will be understood to have an upper summation limit of N, while the index r on x r 
will be understood to have an upper summation limit of 2N.) If for convenience we order 
the bosonic variables before the fermionic ones in the 2N dimensional phase space vector x r , 
then the matrix uj rs which appears in Eqs. (9a, b) is given by 

uj = diag(Q B , ...,n B ,fl F , ....,fl F ) , (10a) 

with the 2x2 matrices Qb,f given by 

Ob=(_° 1 I) , n, = -(; J) , (io 6 ) 

and uj obeys 

(oJ 2 ) rs = -e r S rs , u; 4 = l , u sr = -e r u rs = -e s u rs , 

/ J ^rs^rt — ^sr^tr ~ <J st ■ 

r r 

This concludes our review of generalized quantum dynamics; the reader interested in fur- 
ther details, including the Lagrangian formulation and some concrete examples of models 
constructed using the total trace formalism, should consult Refs. [1-4]. 



2. Weyl Ordered Hamiltonians 

Let us now consider a special class of operator Hamiltonians called Weyl ordered 
Hamiltonians, in which the bosonic operators are all totally symmetrized with respect to 
one another and to the fermionic operators, and in which the fermionic operators are to- 
tally antisymmetrized with respect to one another. Clearly, the most general Weyl ordered 
Hamiltonian which is a polynomial in the operator phase space variables {ay} will be a sum 
of terms, which may be of different degrees, each obtained by Weyl ordering a distinct mono- 
mial in the phase space variables. The contribution of all such monomials of degree n may be 
simply represented by a generating function G n constructed as follows. Let ay, r — 1, 2N 
be a set of parameters which are real numbers when e r — 1 and which are real Grassmann 
numbers, which anticommute with each other and with all of the fermionic phase space 
variables, when e r = —1. Then if we form 

G n = g n , g = ^a s x s , (11a) 

s 

the coefficient of each distinct monomial in the parameters ay will be a distinct Weyl ordered 
polynomial of degree n in the phase space variables {x r }. Corresponding to the operator 
generating function G n , we define a total trace functional generating function 

G n = TrG n , (116) 

where we specify the action of Tr on the Grassmann parameters in G n by requiring that each 
Grassmann oy anticommutes with (— 1) F . The part of G n which is even in the Grassmann 
parameters is then a generating function for all nonvanishing total trace functionals that 
correspond to the bosonic Weyl ordered monomials generated by G n . 

Let us now compare the generalized quantum dynamics equations of motion produced 
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by G n for general operators {x r }, with the corresponding Heisenberg picture equations of 
motion produced by G n when the phase space variables {x r } are assumed to obey the 
canonical algebra of complex quantum mechanics. In our compact phase space notation, 
this algebra takes the form 

(12) 

[x r ,i] = , 

where we adopt the convention that if only one of bosonic, it is taken to be the 

operator x r ; alternatively, we can rewrite the first part of Eq. (12) with no restrictions on 
the indices r, s by including a factor a s , giving 

[x r ,a s x s ] = iu rs a s . (13) 

Applying the equations of motion of Eq. (9b) with G n playing the role of the total 
trace Hamiltonian, we get 

x r = ^u rs ^^ = ^2uj rs ng n - 1 a s . (14) 

s s s 

On the other hand, from the canonical algebra of Eq. (13) we find, for both bosonic and 
fermionic x r , that 

[x r ,g] = i^2uj rs a s , (15a) 

s 

which in turn implies that 

[x r , G n ] = ng n ~ x i ^ u rs a s . (156) 



s 



But the Heisenberg picture equations of motion for the phase space variables, taking G n as 
the operator Hamiltonian, are 



which substituting Eq. (15b) becomes 



x r = to rs ng n 1 a s 



(16a) 



s 



in agreement with Eq. (14). We can now sum over all generating function contributions G n 
weighted by c-number coefficients to obtain a general Weyl ordered Hamiltonian H, which 
has a corresponding total trace Hamiltonian H = TrH, which respectively generate the 
Heisenberg picture equation of motion 



and the corresponding generalized quantum dynamics equation of motion of Eq. (9b). 

Thus, for Weyl ordered Hamiltonians formed with c-number coefficients, we conclude 
that the generalized quantum dynamics equations of motion generated by H agree with the 
Heisenberg picture equations of motion generated by H, on an initial time slice on which 
the phase space variables are canonical, and that on this time slice 



But since Eq. (16c) guarantees that the Heisenberg picture equations of motion preserve 
the canonical algebra on the next time slice, integrating forward in time step by step then 
implies that generalized quantum dynamics agrees with Heisenberg picture dynamics at all 
subsequent times, and therefore defines a unitary dynamics. We have given the argument 
here in a form which applies in complex Hilbert space, where a c-number is a general complex 
number; in quaternionic Hilbert space i in the above equations is actually an operator in 
the left quaternion algebra, denoted by / in [1, 2], and the most general c-number is a real 
number; while in (even dimensional) real Hilbert space % is represented by a 2 x 2 real matrix 



(166) 



[H,i] = 



(16c) 
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i 2 , which has the form of — Qb of Eq. (10b), and the most general c-number is again a real 
number. 

When H is not Weyl ordered, one can give explicit examples in which the generalized 
quantum dynamics equations of motion do not agree with those computed from Heisenberg 
picture quantum mechanics. With one pair of bosonic phase space variables q, p, for example, 
and the fifth degree Hamiltonian 

H =7i(j»V + p 2 q 2 p + pq 2 p 2 + q 2 p 3 + qp 3 q) 

+72 (p 2 qpq + pqpqp + qpqp 2 + pqp 2 q + qp 2 qp) , 
explicit calculation gives 

ATT 

[iH, q] - = 2( 7 i - 72) , (176) 

which vanishes only in the Weyl ordered case 71 =72. The nonvanishing contribution to the 
right hand side of Eq. (17b) can be traced to the operator rearrangement 

(q 2 p 2 — qp 2 q) + (p 2 q 2 — qp 2 q) = Hqp — 2ipq = — 2 , (17c) 

which involves two successive applications of the canonical commutator [q,p\ = i. Taking 
the trace of the left and right hand sides of Eq. (17c) clearly leads to a contradiction if cyclic 
invariance of the trace is assumed. This example serves as a warning that, even though 
taking Tr of Eq. (12) does not directly lead to an inconsistency (because while the cyclic 
property of the trace implies that 

) = Tr[x r x s - (e r ) 2 x r x s ] = , (17d) 

the inclusion of the real part Re in the definition of Eq. (1) makes Tri = 0), the canonical 
algebra is in general inconsistent with the cylic trace property. Hence we cannot simply 
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impose the canonical algebra by fiat as an initial condition in generalized quantum dynamics! 
But we shall see that an effective canonical algebra can arise as an emergent property of 
ensemble averages in the statistical mechanics of generalized quantum dynamics. 

To complete the analysis of this section, we note that Weyl ordering is not a necessary 
condition for the two forms of dynamics to agree, as can be seen by considering the fourth 
degree Hamiltonian 

H = 7! (pV + g V) + i2pq 2 p + i3qp 2 q + ii{pqpq + qpqp) , (I7e) 

for which one finds that the two forms of dynamics coincide for arbitrary values of the coef- 
ficients 7i,...,4 which multiply distinct self-adjoint combinations of operators. This example, 
and analogs with more than one degree of freedom, are relevant for the behavior of the oper- 
ator gauge invariant extensions of standard gauge theories formulated in Refs. [1, 2], which 
we will study in detail elsewhere. 

3. The Conserved Operator C 

Let us now make a further application of the generating function G n of Eq. (lib). 
Taking the operator derivative with respect to x s and using the fact that a s commutes with 
g, we have 

n -i 



(18) 

n-l 



Sx s 

=a s ng 

Multiplying the first equality in Eq. (18) by x s from the right and summing over s, we get 

^2 ~^r x s = n 9 n ~ l ^2 asXs = ngn ' ( 19a ) 

s s s 

while multiplying the second equality in Eq. (18) by e s x s from the left and summing over s, 
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we get 



^2^x s ^- = ( ^2e s x s a s ) ng n 1 = ( ^a s x s \ng n 1 = ng n . (196) 



5x s 

s \ s / \ s 



Since the right-hand sides of Eqs. (19a, b) are the same, subtracting them gives 



s 



which when summed with c-number coefficients over all monomial contributions to the Weyl 
ordered total trace Hamiltonian H yields the important identity 

^{■^'-^'b.) = ■ (20) 

The identity of Eq. (20) is in fact more general than is suggested by the preceding 
derivation, and holds even if H = TrH is not Weyl ordered, provided only that H is con- 
structed from monomials formed from the {x r } using only coefficients that commute with 
all bosonic operators in Hilbert space (that is, as discussed above, real coefficients in quater- 
nionic and real Hilbert space and complex coefficients in complex Hilbert space; in addition 
to these c-numbers, the coefficients can also depend on the grading operator (— 1) F .) To see 
this, let us consider two distinct variations of H which we label Si and <5 2 , defined respectively 
by 

S\x r =x r SA , 

(21) 

^2^r =(SA)x r , 

with 5A an arbitrary infinitesimal anti-self-adjoint bosonic operator variation. As long as 
the only noncommutativity in H arises from the phase space variables {x r }, which is the 
case when only c-number coefficients are employed in constructing H, cyclic invariance of 
Tr implies that the two variations give the same result when applied to H, since the term 
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OlSAOr arising from 5± acting on the right-most operator in Ol is identical to the term 
arising from o~ 2 acting on the left-most operator in Or. Applying Eq. (4) to Eq. (21) gives 

0=Tr]T— [x s 6A-(5A)x a ] 

5H 



<m 

Sx s 

m 

Sx 



Xs ^-sXs 

S VJj s 



' Sx, 



(22a) 



SA , 



which since SA is an arbitrary anti-self-adjoint bosonic operator implies the anti-self-adjoint 
operator relation 



<5H SH 

~r Xg € s Xg — 

OX. Sx, 



(226) 



giving the same identity as was obtained in the Weyl ordered case in Eq. (20). 

We shall now rewrite the identity of Eq. (22b) in an alternative useful form. Let us 
define the operator C by 

^^^^^ OCfUJfgOb g 

(23) 

r,B r,F 

that is, C is the difference between the sums of bosonic commutators and fermionic anticom- 
mutators. (The tilde is a reminder that C is anti-self-adjoint, a point that will be discussed 
in detail shortly.) Differentiating the first line of Eq. (23) with respect to time, and using 
the Hamilton equations of Eq. (9b) and the properties of u) rs summarized in Eq. (10c), we 
find 



(-j — ^^^^^ (^X f LUy g 0C g | 0C i-pUJ i-p gOC g ^ 



r,5 



£ 

r,s,t 

£ 



SH 

Sx L 

SB. 

Sx, 



■Ul r tUJ rs X s ~\- X r 0J rs 0J 



SB\ 



st 



- 



5x t J 



(24) 



SB\ 



- 



5x r J 
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where in the final equality we have used the identity of Eq. (22b). Thus the operator C 
is a constant of the motion for generalized quantum dynamics, as long as the total trace 
Hamiltonian H is constructed from the operator phase space variables using only coefficients 
that commute with all bosonic phase space operators. 

We will exploit the conservation of C in the following two sections, but pause to 
make two remarks. 

(1) We first comment on the adjointness properties of C. Taking the bosonic coordinates 
{q r } and conjugate momenta {p r } to be self-adjoint operators, the bosonic commutator 
terms in Eq. (23) defining C are evidently anti-self-adjoint. (In Refs. [1, 2] the possibility 
of anti-self-adjoint bosonic {q r } and {p r } was considered, and this adjointness assignment 
also leads to an anti-self-adjoint C; however, because of the complex quantum mechanical 
structure of the final results of this paper we expect the anti-self-adjoint case to resemble 
the self-adjoint case, and do not consider it further.) In complex Hilbert space, where i is 
a c-number, the usual fermionic Lagrangians lead to the identifications (see Sec. 13.6 of 
Ref. [2]) q r = ifj r , p r = iipl, with tp r a fermionic operator which is neither self-adjoint nor 
anti-self-adjoint. This gives gj = — ip r , pi = —iq r , as a consequence of which the fermionic 
anticommutator terms in Eq. (23) are also anti-self-adjoint. An analogous construction 
also applies in quaternionic Hilbert space, with the left algebra operator / playing the role 
of i. One must now pay attention to factor ordering and use of a manifestly self-adjoint 
Lagrangian, giving for each fermionic degree of freedom a pair of phase space operators 
q rl = ip r , p r i = q r -2 = and p r2 = fy, so that p\ 2 = q rl , q\ 2 = -p rl , which 

makes {q r i, p r i} + {<2V2, Prt} anti-self-adjoint. In quaternionic Hilbert space, the only way to 
construct a total trace Lagrangian for fermions without using an explicit imaginary unit is 
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to introduce the fermions in pairs, with the real matrix —fig introduced above playing the 
role of the imaginary unit (see Sec. 13.7 of Ref. [2]). For each r one then has q r i,p r i, <?r2, Pr2, 
with adjointness properties assigned according to p\ 2 = ^Qri , q\-2 = TPri (but now with 
no relation between p ri and q ri ), which again makes the two-term sum {q r i,p r i} + {q r 2,Pr2} 
anti-self-adjoint. This fermionic construction also applies to real Hilbert space. 
(2) Not surprisingly the conservation of C can be given a Noether's theorem formulation in 
terms of the total trace Lagrangian L = TrL[{q r }, {q r }] which corresponds to H. We discuss 
here the simplest case, in which L is constructed from its operator arguments using only 
c-number coefficients, and involves no constraints. Cyclic invariance of Tr then implies that 
L is invariant under the operator variations 5q r = 5 2 q r — Siq r = [<5A, q r ] , r — 1, N for 
a time-independent bosonic variation 5 A. The generalization of Noether's theorem to total 
trace Lagrangians given in Sec. 13.5 of Ref. [2] then implies that there is a conserved charge 
Q\ obeying Q A = and given by 



which on substituting p r = 5~L/5q r becomes identical to C. If the label r is a composite 
label comprising a spatial coordinate x as well as a discrete field index r, then the Noether's 
theorem argument implies that there is a current J M which obeys <9 M J M = and is given by 




(25) 




(26) 



with C = J d 3 x J° the associated charge. 
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4. Canonical Transformations 

We turn now to an analysis of the structure of canonical transformations and sym- 
metry transformations in generalized quantum dynamics. Generalizing from the structure [5] 
of infinitesimal canonical transformations in classical mechanics, an infinitesimal canonical 
transformation in generalized quantum dynamics is defined by 

. ^ — > 5G 

x r — x r = ox r = > uj rs - — . (27) 

OXg 

s 

Here G = TrG, with G self-adjoint, is a total trace functional constructed from the phase 

space operators {x r } using arbitrary coefficients, which can be fixed operators as well as 

c-numbers. When Eq. (27) is restricted by requiring that the coefficients used to form G are 

composed of either c-numbers, Grassmann c-numbers, or the grading operator (— 1) F , the 

transformation will be termed an intrinsic canonical transformation, and when the further 

condition of a Weyl ordered G is imposed, the transformation will be termed a Weyl ordered 

intrinsic canonical transformation. 

Letting A = A[{x r }] be an arbitrary total trace functional, we find immediately 

that to first order under a canonical transformation, 

A' =A[K}] 

=A + Tr V^&r r 
^ 5x r 

(28a) 

. m 5 A 5G 

r,s 

=A + {A,G} , 

that is, 

5 A = A' - A = {A, G} . (286) 

Comparing Eq. (27) with Eq. (9b), we see that when G is taken as Hdt, with H the total 
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trace Hamiltonian and dt an infinitesimal time step, then 5x r = x r dt gives the small change 
in x r resulting from the dynamics of the system over that time step. From Eq. (28a), we 
see that when A = TrA with A a Weyl ordered polynomial in the arguments {x r }, the 
canonically transformed total trace functional A + {A, G} obtained by applying a canonical 
transformation is again a Weyl ordered polynomial in the new arguments {x' r }. However, 
since we shall see in Appendix F (where we discuss further details of canonical transforma- 
tions) that the class of Weyl ordered total trace functionals is not closed under the generalized 
Poisson bracket operation, the transformed functional A' is not in general a Weyl ordered 
polynomial in the original arguments {x r }. 

In Sec. 2 we saw that for Weyl ordered Hamiltonians, there is a close relationship 
between the time evolution under generalized quantum dynamics and the corresponding 
Heisenberg dynamics generated when the operator variables are assumed to obey canonical 
commutators. Generalizing the calculation of Eqs. (14-16b) to the case when H is replaced 
by any Weyl ordered intrinsic canonical generator G [the use now of Grassmann c-number 
coefficients in adding monomials causes no problems since g in Eq. (11a) is bosonic], we see 
that in this case Eq. (27) can be represented over the canonical algebra by a commutator as 
well as by an operator derivative, 



with the first equality in Eq. (29) holding for arbitrary operator arguments, and the second 
equality holding only over the canonical algebra. 

In the remainder of this section we establish two important invariances under canon- 
ical transformations in generalized quantum dynamics. We consider first the change in the 




s 



(29) 



s 
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conserved operator C under an intrinsic canonical transformation, giving by use of Eq. (27) 

8C = ^^[5x r uj rs x s + x r u rs 5x s ] 

r,s 

~ 1^ ( T^T^rt^rsXs + X r U rs U st — I (39) 



- V ( — ^ — 

> I — X r G r X r — 

\ ox r ox r 

r x 

where we have again used the identities of Eq. (10c). We recognize that the right hand side of 
Eq. (30) has the same structure as was encountered in Sec. 3, with G now playing the role of 
H. Therefore by the same cyclic invariance argument as was used previously in Eqs. (21-22), 
we conclude that as long as only c-numbers, Grassmann c-numbers, or the grading operator 
(— 1) F , all of which commute with an arbitrary bosonic 5A, are used in constructing G, the 
right hand side of Eq. (30) vanishes and C is intrinsic canonical invariant. 

The second invariance concerns the phase space measure for the phase space opera- 
tors {x r }. Let us introduce a complete set of states {\n)} in the underlying Hilbert space, 
so that the phase space operators are completely characterized by their matrix elements 
(m\x r \n) = (x r ) mn , which have the following form in real, complex, and quaternionic Hilbert 
space: 

In real Hilbert space : (x r ) mn = (x r )° mn 

In complex Hilbert space : (x r ) mn = {x r )° mn + i(x r )^ n (31) 

In quaternionic Hilbert space : (x r ) mn = (x r )° mn + ?(a; r )^ n + j(x r ) 2 mn + k(x r ) 3 mn , 
with (x r )^ n , A = 0, 1,2,3 real numbers. (Note that this is true for fermionic as well as 
bosonic operators; the matrix elements of fermionic operators are still real numbers, not 
real Grassmann numbers! Grassmann numbers are employed only as auxilliary quantities in 
forming Weyl ordered products and in grade-changing transformations.) If for the moment 
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we ignore adjointness restrictions, the phase space measure is defined by 

d/i = J^J dfi A , 



d\x A = ] [ d(x r ) 



A 

mn i 



(32) 



when adjointness restrictions are taken into account, certain factors in Eq. (32) become 
redundant and are omitted. Our strategy is to first ignore adjointness restrictions and to 
prove the canonical invariance of each individual factor dfi A in the first line of Eq. (32); in 
Appendix A we describe how the adjointness restrictions modify Eq. (32) and show that the 
proof remains valid when these modifications are taken into account. 

Under the canonical transformation of Eq. (27), the matrix elements of the new 
variables x' r are related to those of the original variables x r by 

WL = (*r)L + Y, U r.(¥P\ A ■ (33) 
s \ 0X s/ mn 

Inserting a complete set of intermediate states into the fundamental definition 

5G = TrJ2 t — Sx s , (34) 



Sx s 

s s 



we get 

5G = ^C^rY V x ')™ > (35) 

s,m,n,A \ 0X s/mn 

where e m = 1(— 1) according to whether the state \m) is bosonic (fermionic), and where 
e° = 1 and e A = —1, A = 1, 2, 3. (In Refs. [1, 2] the factor of e A was inadvertently omitted, 
but this does not affect the alternative proof of the Jacobi identity given there.) Thus, we 
see that 

fSG\ A _ A dG 

UJ™ m d(x s )£ m ' (36) 
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allowing us to rewrite Eq. (33) in terms of ordinary partial derivatives of the total trace 
functional G, 

<9G 

(K)L = ML + $>rse m e A (37) 

s u \ x s) n m 

Differentiating Eq. (37) with respect to (x r /)f n , n ,, we get for the transformation matrix 

d(, X 'r)mn _ x r r , A ^ 2 G /oo\ 

dMi, n , ~ d -' d -' d -' + 2. ^ 0(*.)^(*rO^ • (38) 

Since for an infinitesimal matrix 5X we have det(l + 5X) ^ 1 + TrSX, we learn from Eq. (38) 

that the Jacobian of the transformation is 

J =1 + 2 , 



y _ d*G • (39) 

r,s,m,n U \ X s)nm U K x r )mn 

Interchanging in E the summation indices r and s, and also interchanging the summation 
indices m and n, we get 

r,«,m,n 1 r ^n^\ x s ) nm 

but now using oj sr = —e r uj rs together with the relation e r = e m e n [which expresses the fact 
that bosonic (fermionic) operators can only connect states of like (unlike) fermion number], 
we obtain 

S = - J2 ^ tmtA d{x )Ad (x )A ■ ( 41 ) 
r,s,m,n °\ X r ) mn O{X s ) nm 

But since the order of second partial derivatives with respect to real matrix elements is 
immaterial, this is just the statement E = — E; hence E vanishes and the Jacobian of the 
transformation is unity. Although we have ignored adjointness restrictions in this argument, 
as shown in Appendix A the conclusion is unaltered when these are taken into account. 
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To summarize, we have shown that the operator phase space integration measure d\i 
is invariant under canonical transformations. An important corollary of this result follows 
when G is taken as the generator dtH of an infinitesimal time translation, since we then 
learn that d\i is invariant under the dynamical evolution of the system, giving a generalized 
quantum dynamics analog of Liouville's theorem of classical mechanics. Since no restrictions 
on the form of the generator G were needed in the above argument for the invariance of 
dfj,, the argument applies even when G is formed from the operator phase space variables 
using operator coefficients. Thus, the integration measure d\i is invariant under a unitary 
transformation on the basis of states in Hilbert space, the effect of which on the variables {x r } 
can be represented by Eq. (27) with G = — Tr J2 r [G,p r ]q r , with G a fixed bosonic anti-self- 
adjoint operator. This transformation, however, is not an intrinsic canonical transformation 
and is only a covariance, rather than an invariance, of the operator C. 

5. Equilibrium Ensemble of Operator Initial Values 

The operator equations of motion of generalized quantum dynamics determine the 
time evolution of the operator coordinates and momenta at all times, given their values on 
an initial time slice. However, these initial values are themselves not determined. We shall 
now make the assumption that for a large enough system, the statistical distribution of initial 
values can be treated by the methods of statistical mechanics. Specifically, we shall assume 
that the a priori distribution of inital values is uniform over the operator phase space, 
so that the equilibrium distribution is determined solely by maximizing the combinatoric 
probability subject to the constraints imposed by the generic conservation laws. Liouville's 
theorem implies that if the assumption of a uniform a priori probability distribution is made 
at one time, then it is valid at all later times, assuring the consistency of the concept of an 
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equilibrium ensemble. We do not propose to address the question of how the randomness in 
the initial value distribution arises: It could come from a random initial condition, an ordered 
initial condition followed by evolution under an ergodic dynamics, or some combination of 
the two. 

More specifically, let dp = dp[{x r }] denote the operator phase space measure dis- 
cussed in detail in the preceding section. In what follows we shall not need the specific 
form of this measure, but only the properties that it obeys Liouville's theorem, and that the 
measure is invariant under infinitesimal operator shifts 5x r , that is 



(This property will not be used until Sec. 6, where we discuss the equipartition or Ward 
identities.) For a system in statistical equilibrium, there is an equilibrium distribution of 
operator initial values p[{x r }], such that 



is the infinitesimal probability of finding the system in the operator phase space volume 
element dp, with the total probability equal to unity, 



The first task in a statistical mechanical analysis is to determine the equilibrium distribution 
P- 

Since equilibrium implies that p — 0, the equilibrium distribution can only depend 
on conserved operators and total trace functionals. In the generic case for a Lorentz invariant 
system, the only conserved operator is C and the only conserved total trace functionals are 



dp[{x r + 5x r }} = dp[{x r }} 



(42) 



dP = dp[{x r }]p[{x r }} 





(436) 
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the total trace Hamiltonian H = p°, the total trace three momentum p, and the total trace 
angular momentum J. However, because the graded trace functionals are all indefinite in 
sign, standard statistical methods lead to a divergent partition function in the generic case. 
We shall therefore restrict our discussion to total trace Hamiltonians H = TrH for which 
the generalized quantum dynamics equations of motion imply conservation of the ungraded 
trace functional H = TrH = ReTrH as well as conservation of H; we shall see in Sec. 7 
below that a large class of models has this property. These models are characterized (see 
Appendices C and G) by having an additional conserved operator F, which when restricted 
to the canonical algebra corresponds to the operator for the conserved fermion number F . 
In addition to the conservation of H, we shall assume the more restrictive condition that 
H is bounded from below; conditions for achieving this are also discussed in Sec. 7. We 
shall also assume henceforth an ensemble which is translation invariant, rotation invariant, 
and Lorentz invariant. Since C is invariant under intrinsic canonical transformations, it is 
Lorentz invariant, and so the equilibrium distribution depends on C; similarly, since F is 
obtained from C by a reordering of fermion factors it is Lorentz invariant as well, and so 
the equilibrium distribution also depends on F. In the ensemble rest frame the equilibrium 
distribution can also depend on H and H (since these are the ungraded and graded mass 
functionals in the rest frame), giving the general equilibrium distribution 

p = F, H, H) . (44a) 

The analysis of the implications of this general equilibrium distribution entails considerable 
algebraic complexity, the full details of which will be published elsewhere. However, as 
summarized in Appendices C and G, the results are all qualitatively similar to those obtained 
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from the simpler equilibrium distribution 

p = p((7,H,H) , (446) 

in which the F dependence is dropped. We shall focus on this simplified case in the exposition 
that follows. 

In addition to its dependence on the dynamical variables, p can also depend on 
constant parameter values, with the functional form of p and the values of the parameters 
together defining a statistical ensemble. Including an anti-self-adjoint operator parameter 
A (which corresponds to the structure of C) and real number parameters f and r, which 
correspond to the structure of H and H, the general form of the equilibrium ensemble 
corresponding to Eq. (44b) is 

p = p(C,A;H,f;H,r) . (44c) 

In the canonical ensemble, we shall see that the dependence on C and A is only through the 
single real number TrXC, and so specializing to this case, Eq. (44c) becomes 

p = p(Tr\C; H, f; H, r) . (44d) 

We shall now show that some significant consequences follow from the general form of 
Eq. (44d), together with the fact that the real function p on the right hand side of Eq. (44d) 
is constructed from its real number arguments using only real number coefficients, and the 
assumption that H and H are constructed from the operators {x r } using only c-number 
coefficients and the grading operator (— 1) F . For a general operator O, let us define the 
ensemble average (0)av by 

(0 )jt y = i*?° . (45a) 
J dpp 
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Then when O is constructed from the {x r } using only c-numbers and (— 1) F as coefficients, 
the ensemble average {0)av must have the form 

{0) AV = F (X,(-1) f ) , (456) 

with the function Fq constructed from its arguments using only c-number coefficients (in 
which we include the f and r dependence). This further implies that both the grading 
operator (— 1) F and the ensemble parameter A commute with (0)av, 

[\,(O) AV } = [(-l) F ,(O) AV } = . (45c) 

Let us now exploit the fact that the anti-self-adjoint operator A can always be 
diagonalized (or, in the case of an even dimensional real Hilbert space, reduced to 2 x 2 
diagonal blocks), by a unitary transformation on the basis of states in Hilbert space, which 
we have seen is also an invariance of the integration measure d/i. The functional relationship 
between A and (C)av then implies that (C)av is diagonal (or block diagonal) in this basis as 
well. As described more fully in Appendix B, this brings {C)av into the following canonical 
form in real (when even dimensional), complex, and quaternionic Hilbert space, 
(C)av = ieffD , i e ff = -i\ ff , i 2 eff = -1 , 

(46a) 

[i e ff,D] = , D real diagonal and nonnegative . 
Although the case of general D is interesting, we shall restrict ourselves in this paper to 

the special case in which D is a real constant times the unit operator; in other words, we 

are assuming that the ensemble does not favor any state in Hilbert space over any other. 

Benefiting from some prescience, we denote this real constant by H, and so have 

(C)av =i e fffr ; 

(466) 

{i eff , (C) AV } =~2h. 
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We turn now to the calculation of the functional form of p in the canonical ensemble, 
which is the ensemble relevant for describing the behavior of a large system which is a 
subsystem of a still larger system. The form of p is determined [6, 7] by minimizing the 
negative of the entropy, 



-S 



I 



dpplogp , 



(47a) 



subject to the constraints 



(476) 



J dpp =1 , 
J dppC =(C)av , 
dppfrH ={TrH) AV , 
J dppTrH =(TrH) AV . 

The standard procedure is to impose the constraints with Lagrange multipliers by writing 
jF=y dpplogp + 9 J dpp + J dppTrXC + r J dppfrH + r J dppTrH , (48a 



and minimizing J 7 , treating all variations of p as independent. (Note that it makes no 
difference whether the constraint for C is introduced through a graded or an ungraded trace; 
the difference is a factor of (— 1) F which can be absorbed into the definition of A.) In order 
for T to have a minimum, it must be bounded below; we shall assume that this is the case 
for sufficiently large r (at a minimum, one needs f > \t\, so that the coefficients of the trace 
of H over the bosonic and fermionic subspaces, proportional respectively to f + r and f — r, 
are both positive). Varying Eq. (48a) with respect to p then gives 



p = exp(-l -9- TrXC - rTrH - rTrH) , 



(486) 
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which on imposing the condition that p be normalized to unity gives finally 

p =Z- X exp(-TrXC - rTrH - rTrH) , 

f „_ . (48c) 

Z= dpexp(-TrXC - rTrH - rTrH) . 

From Eq. (48c) we can easily derive some elementary statistical properties of the 

equilibrium ensemble. For the entropy S, we find 

S = ~ J dpp\ogp = \ogZ + Tr\(C) A v + r(frH) AV + T(TrH) AV . (49a) 

Since the ensemble averages which appear in Eq. (49a) are given by 

5\ogZ 



{C)av = 



sx 



(TrH) AV = 



df 
dlogZ 



dr ' 
Eq. (49a) takes the form 

S^Z-^J^-^-r^l. (49c) 
SX or Or 

Thus the entropy is a thermodynamic quantity determined solely by the partition function. 

Taking second derivatives of the partition function, we can derive thermodynamic formulas 

for the averaged mean square fluctuations of the conserved quantities C, H = TrH, and 

H = TrH, 



Al pd =((TrPC - (TrPC) AV ) 2 ) AV = ({TrPCf) AV - (TvPC) 2 AV = (Trpi) 2 log Z , 

oX 

A% =((H - <*W>av = (&) AV - (H)\ v = ° 2 1<)gZ 



=((H - (H) Av f) AV = (H 2 ) AV - (K) 2 AV = 



(dr) 2 ' 
d 2 logZ 

(49d) 



with P an arbitrary fixed anti-self-adjoint operator. Equations (49a-d) show that the en- 
tropy, the expectations of C, H, and H, and the mean square fluctuations of the latter three 
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quantities, are all extensive quantities which grow linearly with the size N of the system. 



H vanishes as N l l 2 in the limit iV — > oo, and justifies using mean values in imposing the 
constraints in Eq. (48a). 

In the Ward identity derivation of the following section, the distribution p enters 
under a phase space integral in two ways. There are a number of terms in which p appears 
simply as a weighting factor; these terms appear to be dominant and in them we assume 
that the distribution is sharp enough so that unvaried factors of the conserved extensive 
quantity C can be replaced by the ensemble average (C)av, an approximation that should 
become exact in the N — > oo limit. In addition there is a term involving the variation Sp of 
the equilibrium distribution, which gives corrections to the Ward identity that we presume 
to come from very high energy physics. The correction term is evaluated using the formula 
(the normalization factor Z is not varied in the following equations, since we are interested 
here only in variations for which SZ = 0), 



This implies that the ratio of the root mean square fluctuation to the mean for C, H, and 



5p =p5 log p 




(50a) 



and from Eq. (48c) we find for the variations of log p, 



Slogp 

5C 
dlogp 

dH 

dlogp 
dH 



T 



T 



A, 



(506) 



The variations 5 Xs C, S Xs H, and 5 Xs lI corresponding to the variation 5x s are computed in 
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Appendix C, with the results 

^ X ^^^^^ (jJ rp £ (yX j> $X g X g(z<y-X j 

r 

5 Xs H =5 Xs TrH = Tr(-1) F ^ x r Lu rs 5x s , ( 50c ) 

r 

5 Xs lI =5 x TrH = Tr ^ x r uj rs 5x s , 

r 

with Cj rs and a ur = J2 s u) us Cj rs given in Appendix C. In evaluating the correction term, 
we assume that unvaried factors of the conserved extensive quantities H and H can also 
be replaced by their corresponding ensemble averages, again an approximation that should 
become exact in the N — > oo limit. 

As a final remark, in the derivation of the next section we shall follow the conventional 
practice of introducing, for each phase space operator, an operator source which can be varied 
and which is then set to zero after all variations have been performed. It is convenient to 
define the sources p r so that they are all bosonic and self-adjoint. To couple such sources to 
the phase space operators, we employ the auxilliary real or Grassmann real parameters a r 
introduced in Sec. 2, and couple the source term as Trp r a r x r . When r is a bosonic index, 
the source p r takes a distinct value for each r, since a r x r is already self-adjoint in this case. 
When r is a fermionic index in complex quantum mechanics, the source p r takes a distinct 
value only for each pair of fermionic phase space operators q r , p r , so that each distinct p r 
multiplies the combination <J qr q r + cr pr p r , which (remembering that ql = —ip r ) is self-adjoint 
when the Grassmann parameters are taken to obey a\ r = —ia Pr . When r is a fermionic 
index in real or quaternionic quantum mechanics, the source p r takes a distinct value only 
for each quartet of fermionic phase space operators q r ±, p r ±, q r 2, p r 2, so that each distinct 
p r multiplies the combination <Jq rl q r i + c Prl Pri + ^2^2 + Cp r2 Pr2, which (remembering that 
q\ 2 = TPri, PI2 — =t<?ri) is self-adjoint when the Grassmann parameters are taken to obey 
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<rj r2 = ±a Prl , a^ Pr2 = T&q r i- We shall follow the practice of writing the source term as 
Tr J2 r Pri^rXr), with the parentheses a reminder of this implicit grouping. With the sources 
included, the equilibrium distribution and partition function take the form 

p =Z~ l exp[-Tr ^ p r (a r x r )] exp(-TrAC - f H - rH) , 

, T (51a) 

Z= rf/iexp[-Tr^p r (a r a; r )]exp(-TrAC-fH-rH) . 

J r 

Continuing to use the expression (0)av to denote the average of a general operator over 
the equilibrium distribution of Eq. (51a) which includes sources, the variations of logZ with 
respect to its source arguments are related to the averages of the x r by 



{{a r x r )) A v = - ^° gZ . (516) 
dp r 



6. Ward Identities, Unitarity, and the Canonical Algebra 

In the previous sections we have seen that in generalized quantum dynamics there 
is a conserved operator C, given by the sum of commutators for all of the bosonic degrees 
of freedom minus the corresponding sum of fermionic anticommutators, and that this oper- 
ator plays a role in equilibrium statistical mechanics closely analogous to that played by the 
summed energy of independent degrees of freedom in classical statistical physics. This natu- 
rally suggests the idea that the canonical commutation relations of quantum mechanics may 
arise from a generalized quantum dynamics analog of the classical theorem of equipartition 
of energy. To pursue this thought, let us begin by reviewing a simple derivation [8] of the 
classical equipartition theorem. Let H({x r }) be the classical Hamiltonian as a function of 
classical phase space variables {x r }, and let dp({x r }) be the classical phase space integration 
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measure. We consider the integral 




^ d[x r exp(— 13 H)] 
dx s 

dfj,5 rs exp(— f3H) 
dfix r ^^^ exp(—j3H) 



(52a) 



the left hand side of which is the integral of a total derivative and vanishes when the integrand 
is sufficiently rapidly vanishing at infinity. Assuming this, we get 



which is the classical theorem of equipartition of energy. The method of derivation is similar 
to that used to derive Ward identities from functional integrals in quantum field theory (see, 
e.g. [9]), and the equipartition theorem can be viewed as a Ward identity application in 
classical statistical mechanics. 

We proceed now to derive a Ward identity for the statistical ensemble of generalized 
quantum dynamics. The derivation is based on Eq. (42), which asserts the invariance of the 
operator phase space measure dfi under finite operator shifts 5x r , which can be arbitrary 
apart from the obvious restriction that they must satisfy the same self-adjointness restric- 
tions as the corresponding operators x r . We take account of the adjointness restrictions on 
the variations by using the following Lemma, proved in Appendix D, which insures that when 
we equate the variation of a total trace functional to zero, we do not inadvertently "deduce" 
an operator relation which arises from the variation of an anti-self-adjoint operator, which 
is identically zero when acted on by Tr. 



Let Y\ and Y 2 be two self-adjoint bosonic or two anti-self-adjoint bosonic operators con- 



/ dnx r (3{dH/dx s ) exp(-/3#) 
f dfi exp(—j3H) 



(526) 



Lemma: 
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structed from the phase space variables. Then in = 5TrYiY 2 , the self-adjointness restric- 
tions on the variations can be ignored. 



To derive the Ward identity, we consider the expression 



= / d/j,5, 



(53a) 



expf-TrAC - f H - rH - Tr Pr(<*rX r )] Tr{C, i eff }V 

where the operator variation 5 Xs is defined to act on an arbitrary operator X[{x t }] as 

5 Xs X[{x t }]=X[{x t ,t^s;x s + 5x s }]-X[{x t }] , (536) 

and where in equating the shift to zero we are using the shift invariance of the measure 
and the assumption that the integrals are sufficiently convergent that contributions from 
infinity can be ignored. In Eq. (53a), the expression V denotes any self-adjoint polynomial 
in the variables {(cr s x s )} constructed using coefficients which are c-numbers apart from a 
possible dependence on the operators (— 1) F and i e ff. [Inclusion of the auxilliary factors a s 
in the combination (a s x s ), which was defined in the discussion preceding Eq. (51a), is purely 
a matter of convenience; it permits working with bosonic quantities throughout, and also 
facilitates comparison with the form of the canonical algebra given in Eq. (13).] The traces 
in the exponent in Eq. (53a) and the trace involving V both have the form specified by the 
Lemma, so we can proceed with taking variations, giving [c.f. Eqs. (50a, b)] 
0= / a>exp[-TrA<5-fH-rH-Tr^p r ((X r :r r )] 

J r 

x [-Tr\5 Xs C-r5 x H-r5 x H-Trp s a s 5x s }Tr{C,t eff }V (53c) 

+ Tr({5 Xa C : i eff }V + {C,i eff }5 X3 V) 
We now make two assumptions: First, when the extensive quantity C appears in 
unvaried form as a factor in an ensemble average over the equilibrium distribution, we assume 
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that it can be replaced by its ensemble average {C)av (but this is not, of course, applied to the 
C in the exponent of the equilibrium distribution). This assumption amounts to neglecting 
the fluctuations of C when it appears as a factor in the integrand, and as we argued in 
Sec. 5, should be justified in the limit N oo. Second, we assume the form of Eq. (46b) for 
(C)av- Replacing unvaried factors of C by their ensemble averages, substituting Eqs. (46b) 
and (50c), and making some cyclic permutations under the graded trace, Eq. (53c) becomes 

= J dp expf-TrAG* - f H - rH - Tr p r (o- r x r )] 



(— Tr[A, u rs x r ]5x s — f Tr(— 1) F x r u rs 5x s 
r r 

— rTr x r u rs 8x s — Trp s a s 5x^j (—2h)*V 



(54a) 



5V 



+ Tr[{i e ff, V}, uj rs x r \8x s — 2uTr- — 5x 

where in the final line we have used the definition of the operator derivative of the total trace 
functional V = TW, 

SV 

Tr5 Xs V = Tr—5x s . (546) 

Before proceeding further, let us examine the structure of the first term on the right 
hand side of Eq. (54a), which contains the factor [A, ^ r u rs x r ]. After varying with respect 
to the sources and setting the sources to zero, this term is proportional to 

with p the zero source equilibrium distribution of Eq. (48c) and with P the polynomial 
which results after variation with respect to the sources. Since V, by assumption, involves 
no operator coefficients other than (— 1) F and ? e //, the reasoning of Eqs. (45a-c) implies 
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that the ensemble average in Eq. (54c) is a function G(\, (— 1) F ), and so commutes with A. 
Hence the first term on the right hand side of Eq. (54a) vanishes. Since each remaining term 
in Eq. (54a) is the graded trace of an operator times the variation 5x s , we can equate the 
total operator coefficient of this variation to zero. After multiplying through by \ ^ s uj us , 
and using Eq. (10c) and the definition of a ur following Eq. (50c), this gives the operator 
Ward (or equipartition) identity 

= / dpexp[-Tr\C - tH - rH - Tr^ p r (a r x r )} 

J r 

[f(-l) F ^2a ur x r + tx u + ^2uj us a s p s }h~V (55 a ) 



x 



r l r . , .v- 5V 

+ [TjVeff, V\, X u \-h2_^ UJ US - - 



s s 

Dividing by the partition function Z, Eq. (55a) can be rewritten in the compact form 

1 5~V 

= ([T(-l) F ^a ur x r + rx u + ^LJ us a s p s ]hV+[-{i eff ,V},x u }-h^u us — ) AV , (556) 

r s ~ s 

with ( )av denoting the ensemble average with sources present, and with the understanding 
that after variation with respect to the sources, the sources are to be set equal to zero. A 
detailed discussion of how symmetrized polynomials in the phase space operators may be 
built up through source variation is given in Appendix E, and a discussion of a second 
Ward identity connected to the conserved operator F is given in Appendix G. Although the 
explicit source term in Eq. (55b) proportional to J2 S ^usVsPs does contribute when varied 
with respect to a source p v for which uj uv is nonzero, the fact that cu uv is antisymmetric 
in bosonic indices, and symmetric in fermionic indices, implies that this term drops out of 
Weyl ordered expressions in x u and the additional factors of x v , ... brought down by source 
variation, as is discussed in more detail in Appendix E. Hence we shall drop this term 
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in the Ward identity applications that follow. A closely related remark is that although 
Eq. (55b) makes a statement of effective equality (when the f and r terms are dropped) 
between the commutator expression in the next to last term and the total trace derivative 
in the final term, this does not contradict the counterexample of Eqs. (17a, b) above, since 
the Ward identity does not imply that the commutator can be evaluated in terms of the 
canonical algebra by using the Leibnitz product rule. Use of the Leibnitz product rule for 
V is justified in Appendix E only when V can be constructed as a Weyl ordered polynomial 
in the operator phase space variables {x r }. 

We proceed now to give three applications of Eq. (55b). As our first application we 
choose V to be the operator Hamiltonian H, so that V becomes the conserved quantity H. 
Substituting the generalized quantum dynamics equation of motion of Eq. (9b), the Ward 
identity becomes 

= ([f(-l) F ^2 a urX r + rx u ] hB. + [^{i e ff, H}, x u ] - hx u ) A v • (56a) 

r 

We can simplify Eq. (56a) considerably by noting that since H is a conserved extensive quan- 
tity, in the large N limit we can approximate it by its ensemble average (H)av] its coefficient 
in Eq. (56a) is then proportional to the partition function variation ^2 s u us SZ/5x s = 0. So 
dropping the f and r terms, and multiplying Eq. (56a) by —1 we are then left with 

1 

= {hx u - [-{i eff ,H},x u ]) A v , (566) 

where the ensemble used to form the average is understood to still contain nonzero sources. 

It is convenient at this point to recall the properties of i e ff given in Eq. (46a), and 
to make the definition 

Heff = l{ieff,H}, Hl ff = -H eff , (57a) 
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with the tilde indicating that if e // is an anti-self-adjoint operator. We find from Eq. (57a) 
that 

ieffHeff = \{~H + ieffHi ef f) = H eff i eff , (576) 

in other words, i e ff and H e ff commute. We can now write Eq. (56b) as 

= (hx u - [H e ff,x u ]) AV , (57c) 

which is an effective Heisenberg picture equation of motion for x u in anti-self-adjoint gen- 
erator form. Let us now vary with respect to the sources, leading (as described in Appendix 
E) to the replacement of Eq. (57c) by the expression 

= (hP({x r }) - [H eff , P({x r })]) AV , (57d) 

with P({x r }) a Weyl ordered polynomial formed with coefficients which are c-numbers 
(apart from a possible dependence on (— 1) F and i e ff)- In particular, letting P be the 
effective Hamiltonian i? e //, we find that 

±(H eff ) AV = , (57e) 

and so (H e ff) A v, still in the presence of sources, is a constant of the motion. Thus the 
ensemble averages of products of the coordinates have an effective unitary time development, 
involving an anti-self-adjoint effective time independent Hamiltonian. 

This time development, however, cannot immediately be put into the standard 
Heisenberg picture form of Eq. (16b), which involves a self-adjoint Hamiltonian. We now 
show that we can extract from the {x r } a set of new operators {x r e //}, which do obey an ef- 
fective dynamics of the standard Heisenberg form. We begin by introducing the self-adjoint 
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effective Hamiltonian H e ff defined by 

H e ff = -i e ffH e ff = \{ H ~ i e ffHi ef f) , (58a) 
which evidently also commutes with i e ff. In analogy with Eq. (58a), we further define 

X r eff = -(x r - i eff X r i e ff) , (586) 

which obeys 

ieffXr eff = ^effXr + X r i e ff) = X r e ffi e ff , (58c) 

and thus also commutes with i e ff. For any operators x±, x 2 this definition evidently obeys 

(x 1 X 2e ff)eff = (xieffX 2 )eff = X leff X 2 eff , (58d) 

and so taking the effective projection of the equation of motion of Eq. (57c) gives 

= (hx u e ff - [H eff ,x u effX)Av , (59a) 

which by Eq. (58c) can now be written directly in terms of the self-adjoint effective Hamil- 
tonian defined in Eq. (58a), 

= (kx u eff - i eff [H eff , x u e ff\)Av , (596) 

and so has the standard Heisenberg picture form of complex quantum mechanics. In a 
similar fashion, by repeated applications of Eq. (58c) to polynomials of successively one 
higher degree, we can derive a self-adjoint analog of Eq. (57d), 

= (hP({x r eff }) - i eff [H eff , P({x r eff })])Av ■ (59c) 

Before proceeding to further Ward identity applications, we make two remarks. The 
first is that the method of projecting out effective operators which commute with i e ff is 
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simply a complex analog of the method of extracting "formally real" components of operators 
in quaternionic quantum mechanics [2]. The second is that the operators x r and x r e ff always 
differ in the cases of real and quaternionic Hilbert spaces. Even in the case of a complex 
Hilbert space they differ when i e ff ^ i, as shown in Appendix B, while in Appendix G we 
show that the complex case with i e ff = i is excluded. When x r and x r e ff differ, there is a 
"hidden sector" which cannot be attained by acting with arbitrary polynomials formed from 
the effective operators alone. 

We turn now to the second Ward identity application, which we derive by taking V in 
Eq. (55b) to be V = (a t x t ). The parentheses here indicate summation over the one, two, or 
four t values corresponding to a distinct source term p t in the equilibrium distribution, which 
according to the discussion preceding Eq. (51a) makes V self-adjoint. Since the a parameters 
which appear in this sum are linearly independent, we can ignore the parentheses and simply 
substitute the single term V = a t x t into Eq. (55b), corresponding to 5V/5x s = <7 t 5 ts , giving 

= ([r(-l) F ^a ur x r + tx u ] hTrcr t x t + [^{i eff ,a t Xt},x u ] - HuutCTt )av ■ (60a) 

r 

Referring back to Eqs. (58a-d), we see that Eq. (60a) can be rewritten in terms of x t e ff as 
= ([f(-l) F ^cwr r + Tx u ]hTra t x t e ff + [ieff<?tXt eff,x u ] - hu ut cr t ) A v ■ (606) 

r 

Taking the effective projection of Eq. (60b) by using Eq. (58d), and then multiplying through 
by i ef f, we get finally 

= (i e ff[r(-l) F ^2a ur x r eff + Tx u eff ]hTr(T t Xt eff + [x u eff ,cr t x t eff ] -i eff huj ut a t ) AV , 

r 

(60c) 

which coincides in form with the canonical algebra of Eq. (13) (where we had set h — 1) 
when the dynamics dependent terms proportional to f and r are dropped. Since these two 
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terms are proportional to time derivatives x e ff, their neglect should be justified in a regime 
characterized by energies much lower than the energy scale of the underlying dynamics. 
(Also, since these two terms are proportional to the graded trace Tra t x t e ff, they may be 
further suppressed by boson fermion cancellation.) By repeated applications of Eq. (58d) to 
the cases in which V is a successively one degree higher self-adjoint polynomial, one proves 
similarly from Eq. (55b) that after the f and r terms are dropped, 

= ( [x u e// , V (K e// » ] - i eff h "us SV( £ Xr e//>) )av , (60d) 

s 0X s eff 

corresponding in form (when h=l) to the canonical algebra expression comprising the second 
and third terms of Eq. (29). Similarly, as discussed in Appendix E, by varying the sources 
in Eq. (60b) one can justify (when V is Weyl ordered) the evaluation of the commutator in 
Eq. (60d) in terms of the canonical algebra of Eq. (13) using the Leibnitz product rule. 

As a final application of Eq. (55b), we examine its implications for canonical transfor- 
mations. Replacing V now by G, with G the generator of a Weyl ordered intrinsic canonical 
transformation, comparing with Eq. (27) and dividing by Ti, and dropping the f and r terms, 
we see that Eq. (55b) takes the form 

= ( \^h~ l {i e ff,G},x u ] - 5x u ) AV . (61a) 

Defining the anti-self-adjoint generator G e ff by 

Geff=\^eff,G} , (616) 

Eq. (61a) can be rewritten in the form 

= ( 5x u - fr^Geff^Xv^Av , (61c) 
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which after taking the effective projection gives 

= ( 8x u ef f - H'^Gef^Xu e ff])AV , (61d) 

and corresponds in form (when H — 1) to the first and third terms of Eq. (29). Equations 

(61c, d), and their Weyl ordered polynomial generalizations 

=( 5P({x r }) - h- l [G eff , P({x r })]) AV , 

(61e) 

=( 5P({x r e// }) - h-'iGeff, P({ Xr ef f})})AV , 

which are analogous to Eqs. (57d) and (59c), indicate that there is a correspondence between 
Weyl ordered intrinsic canonical transformations in the underlying generalized quantum 
dynamics, and unitary transformations U e ff of the form 

U eff = exp(^ 1 G e// ) (62) 

acting on the variables {x r e ff}- The invariance of i e // under unitary transformation by 
U e ff is an image, in the effective theory, of the invariance of C under intrinsic canonical 
transformations in the underlying generalized quantum dynamics. 

To summarize, we have shown that there is a striking correspondence between the 
structure of the set of ensemble averages calculated in generalized quantum dynamics and 
the structure of canonical quantum field theory. To what specific field theoretic structures do 
these averages correspond? To answer this question we note that in the absence of sources, 
the averages of monomials constructed from the phase space operators all are functions, 
constructed with real number coefficients, solely of the operator i e ff and of the grading op- 
erator (— 1) F . (This statement follows from an argument given at the end of Appendix B.) 
Since the Hamiltonian H is necessarily bosonic, physical processes cannot change the value 
of (— 1) F , and so the (— 1) F = 1 and (— 1) F = —1 sectors of Hilbert space are superselec- 
tion sectors, the states of which do not superimpose. Within each superselection sector, the 
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ensemble averages of monomials are simply linear combinations with real coefficients of the 
two operators f and i e ff\ this implies that they are to be indentified with matrix elements, 
rather than operators, in an effective complex field theory in which i e ff plays the role of the 
imaginary unit. Since for any self-adjoint operator A formed from the phase space operators 
{x r }, reality and positivity of the equilibrium phase space density p imply that (A)av is re& l 
and (A j< A)av is nonnegative, the ensemble averages must correspond to expectation values 
in some pure or mixed state in the effective field theory. But since the ensemble in which 
the averages are formed is Lorentz invariant (recall that in a general Lorentz frame H gets 
replaced by the invariant [— (Tr p^) 2 ]^ , and similarly for H), and since the vacuum is the only 
Lorentz invariant state in quantum field theory, it is then natural to identify the ensemble 
averages with vacuum expectation values in an effective quantum field theory. This leads us 
to conjecture that in the limit of infinitely many degrees of freedom and an infinite dimen- 
sional underlying Hilbert space, the ensemble averages of Weyl ordered polynomials formed 
from the canonical phase space operators in generalized quantum dynamics are isomorphic 
in structure to the vacuum expectation values of the corresponding Weyl ordered polynomials 
formed from the canonical operators in complex quantum field theory. 

An interesting feature that has emerged from our calculations is that an effective 
complex structure results irrespective of whether one starts from an underlying real, complex, 
or quaternionic Hilbert space. In particular, there now seems to be no reason to exclude the 
aesthetically appealing case of real Hilbert space, since we have given a natural and automatic 
mechanism for the complexification which is needed to describe the observed physical world. 
We also note that there is a natural connection between our calculations and matrix models 
containing 2N matrices, which correspond to our 2N phase space operators, acting on an M 
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dimensional space of states in the large M limit (conventionally called the "large TV" limit 
in the matrix model literature). When the term in the equilibrium ensemble containing A 
vanishes, the ensemble is invariant under unitary transformations, in which case the large M 
limit is known to be described by the classical field theory [10, 11, 12] of a so-called "master 
field" . When the A term is nontrivial it breaks the unitary invariance, raising the interesting 
possibility that the master field in this case is a complex quantum field. 

We make next some observations concerning the use of Weyl ordering in our analysis. 
It is clear from the discussion of Sees. 2-4 and Appendix E that the proposed isomorphism 
must fail for certain types of polynomials which are not Weyl ordered. There is evidently 
a subtlety in the non-Weyl ordered case, an understanding of which will require a more 
detailed investigation. Since the operator gauge invariant models discussed in Refs. [1, 2] do 
not have Weyl ordered Hamiltonians, it is important to extend our derivations so that this 
case (which involves fourth and lower degree polynomials in the phase space operators ) is 
explicitly included; this issue will be addressed in a separate publication. 

Finally, we remark that in generalized quantum dynamics the concept of "tempera- 
ture" is not defined, because the operator Hamiltonian H is not a constant of the motion: 
only its graded and (in the class of models studied here) ungraded traces TrH and TrH are 
conserved. However, in the effective complex theory the effective Hamiltonian is conserved, 
making possible further fine grained equilibria governed by H e ff, and thus leading to the 
emergence of the standard thermal ensemble exp(—(3H e ff) parameterized by the temperature 

r 1 . 

7. Conditions for Convergence 

The Ward identities derived in the preceding section are meaningful only when the 
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phase space integrals which appear in them are convergent. As we have already noted, the 
generic conserved quantities C and H are both indefinite in sign, and so when standard 
statistical mechanical methods are applied to them one finds a partition function which is 
divergent. A necessary condition for convergence is that there be at least one constant of 
motion which is bounded below, and an obvious candidate for this is the ungraded trace 
TrH . A simple criterion can be given for constructing Lagrangians, in which the total trace 
equations of motion obtained using the graded trace Lagrangian TrL imply conservation of 
TrH as well as conservation of TrH. Let us consider Lagrangians in which the fermion fields 
appear only through one of the standard bilinears of the form ipj . . .ip r ; then when bosonic 
variables to which the fermions couple are ordered to the outside of the fermion factors, as in 
the vector coupling ^7°7 At ...ip r B^ or the scalar coupling ipj'y ...ip r B, the ungraded trace and 
graded trace equations of motion are the same and imply conservation of TrH. The reason 
is that when the bosons are ordered to the outside of the fermions, only cyclic permutations 
of bosonic factors are required to construct the equations of motion for the bosons and for 
ifj r , and so these are the same irrespective of whether or not a grading factor is included 
inside the trace. The equation of motion for ip\ can then be obtained as the adjoint of one 
of the ipr equations of motion, or if calculated directly by permuting the factor Sip} to the 
right, the same grading factor appears in all terms and so drops out of the ipi equation of 
motion. Closer examination shows that in the models with a conserved TrH, the equations 
of motion also imply conservation of an ungauged fermion number current F^, with J d 3 xF° 
yielding the conserved charge F discussed briefly in Sec. 5 and in Appendices C and G. 

On the other hand, if the bosonic variables to which the fermions couple are ordered 
to the inside of the fermion factors, as in the vector coupling ^^^B^ipr, then ip r must be 
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cyclically permuted to the left in order to obtain the equation of motion, resulting in 
differing signs for the source terms in the graded and ungraded cases. These features of 
the generalized quantum dynamics equations of motion are readily verified by examination 
of the catalog of models given in Sees. 13.6-13.7 of [2]. In generalizing the Lagrangians for 
the various gauge theory and Higgs components of the standard model, and its grand unified 
extensions, so as to give generalized quantum dynamics Lagrangians, one is always free to 
adopt the ordering convention which leads to conservation of TrH. However, one can readily 
construct generalized quantum dynamics models in which TrH is not conserved; an example 
is the maximally operator gauge invariant model constructed in [1, 2] using two fermion 
fields, in which the gauge potential B' is ordered to the outside of the fermion fields and a 
second gauge potential B^ is ordered to the inside, the difference in ordering being precisely 
what distinguishes the action of the two gaugings. In this case, one readily verifies that the 
source term for B^ changes sign as one goes from the graded to the ungraded equations of 
motion, while that for B'^ has the same sign in both cases. 

In addition to the requirement that TrH should be conserved, it is also necessary 
that rTrH should be bounded below and should dominate over the indefinite terms TrAC* 
and rH, in order for the partition function to converge. In general, TrH contains three types 
of terms: bosonic kinetic energy terms, fermionic kinetic energy terms, and potential terms. 
The bosonic kinetic energy terms in standard field theory models are always positive, and 
we believe that in many models one will be able to establish that the potential is bounded 
below. However, because of negative energy states the fermionic kinetic energy terms are 
problematic. For a single Dirac fermion, before any reorderings of field operators are done, 
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the kinetic energy (including the mass term) has the form 



with m the fermion mass. Since Eq. (62a) is the difference of two positive semidefinite terms, 
TriJkin is unbounded below when the operators bp s and dp tS are independent. However, for 
a self-conjugate Majorana fermion field one has dp s = b^ s , with T the operator transpose, 
and so in this case the trace of the second term in Eq. (62a) is 

and cancels the trace of the first term; thus the trace of the kinetic energy term vanishes and 
the negative energy catastrophe is avoided. Hence in order for a field theoretic model to be 
extendable to a generalized quantum dynamics model with a convergent partition function, 
one must be able to rewrite its fermion kinetic energy terms entirely in terms of Majorana 
fermions, which requires that every chiral fermion be accompanied by an opposite chirality 
partner. To complete the discussion of the fermion kinetic energy terms, let us examine the 
behavior of the graded trace of the fermion kinetic energy. For any two operators 0\ and 
O2, irrespective of their grade, we have 

TrOfOj =Re^e m OL„^2 T nm = Re ^ e m O lnm 2mn 

m.n m.n . . 

' (62c) 

m,n 

hence the rearrangement of Eq. (62b) is also valid inside Tr, and so for Majorana fermions 
the graded trace of the kinetic energy term also vanishes. [The difference between Eq. (62c) 
and the cyclic identity of Eq. (2) is that to interchange the untransposed operators one must 
replace e m by e n , and the product e m e n is just the grade of the operators Oi j2 , leading to the 
± sign in Eq. (2).] 
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Although a discussion of sufficient conditions for convergence of the partition function 
is beyond the scope of this paper, let us give one example which suggests that convergence 
can be attained. Let us work in complex Hilbert space, and consider special simple choices 
for the ensemble parameter A. The simplest possibility, A = iR, with R a positive real 
number, is inadmissible because % is a c-number and so TrAC reduces to a multiple of TrC, 
which is identically zero, 

TrC* = Tr x r u rs x s = Tr x s t s uo rs x r = — Tr x s u sr x r = . (63a) 

r,s r,s s,r 

So let us consider the next simplest case, which is A = (—l) F iR, for which 
TrAC =Re iRTiC 



-Re iRTr(£2[q r ,p r ] - ^{q r ,Pr}) 

r,B r,F 

- - Re iRTrJ2Ur,Pr} = OTr ^{g r ,gj} , 



(636) 



r,F r,F 

where in the final line we have assumed that all fermions have been constructed with the 
standard adjointness assignment p r = iq^. Since both terms in the anticommutator on the 
right hand side of Eq. (63b) are positive semidefinite, the term TrAC' is bounded below in this 
case, and does not have to be dominated by the energy term fTrH. So a model in which the 
ungraded trace energy is conserved, bounded below, and becomes positive infinite on all paths 
approaching infinity in operator phase space (which our preliminary investigations suggest 
may be possible for gauge theories) then gives a convergent partition function. Although in 
this special case the A factor in the exponent has decoupled from the bosons, as long as the 
bosons and fermions interact the breaking of unitary transformation invariance implied by 
the presence of A is still felt by the bosons, and the model is satisfactory. For a noninteracting 
theory, of course, the fermionic and bosonic contributions to C are separately conserved; both 
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must then be included in the statistical mechanical equilibrium distribution, and the choice 
A = (— l) F iR would then lead to classical mechanics, rather than quantum mechanics, for 
expectations defined within the bosonic ensemble. 

8. Discussion 

In the foregoing, we have given strong evidence that the canonical quantization 
rules, which are the basis of conventional quantum mechanics and quantum field theory, can 
arise as an emergent property in a generalized quantum operator dynamics when statistical 
mechanical methods are applied to the ensemble of operator initial values. In models where 
the index r is a composite index composed of a spatial coordinate x as well as a discrete 
field index, the emergent canonical algebra implies locality, even though the underlying 
generalized quantum dynamics is highly nonlocal. To our knowledge, this is the first time 
that an embedding of quantum mechanics in a larger structure has been achieved that applies 
to local relativistic quantum field theories. 

An interesting feature of the Ward identity applications of Sec. 6 is that the derivation 
of a unitary effective dynamics requires a less stringent approximation (the replacement of 
H by its ensemble average) than does the derivation of the canonical algebra (where a "low" 
frequency approximation is needed). Thus, there may be a high energy domain where the 
statistical mechanical analysis is still valid and takes the form of a nonlocal but unitary 
complex effective field dynamics, such as a string-type theory. This effective theory would 
still be only a statistical approximation to the fully nonlocal underlying generalized quantum 
dynamics, and its structure may well prove more intractable than that of the underlying 
theory. 

An exciting aspect of our construction is that the principal features of quantum 
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mechanics basically become statements about the geometry of matrices. Thus, if in labeling 
the rows and columns of matrix operators one orders all the bosonic states before all the 
fermionic states, then the distinction between bosonic and fermionic operators is simply that 
between matrices that are block diagonal, and ones that are block skew diagonal, respectively. 
Furthermore, as we have seen in Sec. 6, the complex structure, canonical algebra, and 
unitary dynamics of quantum mechanics are all reflections of the cyclic invariance of the 
trace, which is both the origin of the conserved operator C and the basis for erecting a 
generalized dynamics on noncommutative phase space. These observations suggest that the 
distinction between matter degrees of freedom on the one hand, and gravitational or metric 
degrees of freedom on the other, may be similarly rooted in some simple geometric property 
of generalized quantum dynamics. 

If our conjectured isomorphism can be proved, and if a generalized quantum dynam- 
ics underlies the observed universe, there will be profound implications for some of the vexing 
issues in conventional quantum mechanics. One of these issues is the quantum measurement 
problem. In the underlying generalized dynamics there are no "dice" : the underlying dynam- 
ics is a generalization of classical mechanics to noncommuting phase space operators and is 
deterministic, although not in general unitary. However, the ability to follow this determin- 
istic evolution in detail is lost at the level of the statistical ensemble average, where a unitary 
conventional quantum mechanics emerges. In this picture, a calculation of corrections to the 
ensemble average approximation should permit the resolution of the troubling "paradoxes" of 
quantum measurement theory. A second issue where our picture has important ramifications 
is the problem of infinities in quantum field theory. These divergences arise, fundamentally, 
because of the singularity of the local canonical commutator/anticommutator structure of 
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field theory, which we have argued is an emergent property of ensemble averages in general- 
ized quantum dynamics. The underlying dynamics is nonlocal and nonsingular, and should 
give finite answers to physical calculations. 

There are clearly a number of important questions that must be addressed in future 
work. One of them is to give a detailed justification of the statistical mechanical aspects 
of our calculation, including a complete classification of Hamiltonians for which one can 
prove convergence of the partition function, a proof that a thermodynamic limit exists, and 
a justification of the assumption that each unvaried factor of C (and of other conserved 
extensive quantities) in the Ward identity can be replaced by the corresponding ensemble 
average. The emergence of quantum mechanics from the Ward identity required the neglect 
of the f and r terms, which may be reasonable only in models that have a very large 
ratio of the high mass scale characterizing pre-quantum mechanical physics to the low mass 
scale characterizing quantum physics, and possibly also a high degree of boson-fermion 
symmetry as well. Thus, finding a model in which neglect of the dynamics dependent 
terms can be justified may be tantamount to finding a model that solves the "hierarchy 
problem" of explaining the extraordinarily large ratio between the Planck mass scale and the 
standard model mass scale in particle physics. Since the entire generalized quantum dynamics 
formalism seems to naturally invite the incorporation of boson-fermion symmetries, such as 
generalized forms of supersymmetry, it will be important to analyze such symmetries. A 
possibly related question is to determine what happens when the numbers N B , N F of bosonic 
and fermionic degrees of freedom are equal. This issue arises because if it were evaluated by 
restriction to the canonical algebra, the conserved operator C would be equal (with H—l) to 
i e ff(N B — N F ), which vanishes when N B = N F . This suggests that the conditions for validity 
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of our analysis may be particularly delicate - and interesting - in models with a high degree 
of boson-fermion symmetry. Yet another question is what happens in constrained theories, 
such as the operator gauge invariant theories formulated in Refs. [1, 2]. In canonical gauges, 
where the constraints can be explicitly eliminated, our analysis should apply directly to the 
Hamiltonian restricted to the constraint surface, which involves only the physical degrees 
of freedom. However, there is likely to be an analog in generalized quantum dynamics 
of the methods used to treat constrained quantum field theories (such as Faddeev-Popov 
determinants and BRST invariance), which would permit working in non-canonical gauges 
as well, and it would be interesting, and possibly important, to find it. Beyond these basically 
technical questions is of course the larger issue of whether one can incorporate gravitation 
in a natural way, and whether one can find a compellingly beautiful total trace dynamics 
which gives rise to all of the observed forces and matter fields. 
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Appendix A. Inclusion of Adjointness Restrictions in the 
Argument for Invariance of the Phase Space Measure 

Inspection of the argument of Eqs. (39-41) shows that the diagonal (m = n) and 
off-diagonal (m ^ n) terms in the sum £ vanish separately, and for each of these, the 
summed contribution from the canonical coordinate and momentum pair q r ,p r for each fixed 
r also vanishes separately. This observation permits us to take the adjointness restrictions 
into account; in the following discussion we shall write dfx = d/isd/iF, with d\iB and d\iF 
respectively the bosonic and fermionic integration measures. There are three cases to be 
considered: 

(1) For a bosonic pair of phase space variables q r ,p r , the x r variables are independent but 
are both self-adjoint, and thus 

( X r)mn = eA ( X r)nm ■ (Ala) 

This means that the integration measure must be redefined to include all diagonal terms 
in m, n, but only the upper diagonal off-diagonal terms, so that the bosonic integration 
measure becomes 

dfJ>B = IJ dfig , 

ir ( Alb ) 

d^ = I] d ML ■ 

r,m<n 

The argument for the diagonal terms in this product proceeds just as did that for the di- 
agonal terms in the unrestricted case, while the argument for the off-diagonal terms uses 
Eq. (Ala) in place of an interchange of the summation index pair m, n, together with the fact 
that for a boson e r — 1, again leading to the conclusion that the diagonal and off-diagonal 
contributions to E vanish independently. 
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(2) For a fermionic pair of phase space variables in complex Hilbert space constructed ac- 
cording to the recipe q r — ip r , p r = iipl = iql, the x r variables are no longer independent. 
However, this construction implies that 

(lr)m n = (Pr) nm > (Pr)mn = (lr) nm > (A2) 

and thus, in a complex Hilbert space, the fermionic integration measure must be redefined 

as 



n TT n (A3a) 
dn F = J [ d(x r ) mn 



Similarly, for the analogous fermionic construction p r i = ■■■ in a quaternionic Hilbert 

space, the fermionic integration measure must be redefined as 

d\iF = ] [ dft 

A=0,2 
djlp = J\ d(x r i) 



1>F , 

A -" 3 A (A31) 



Ivan 



Since the argument for the unrestricted case worked for each A value separately, it still goes 
through as before. 

(3) Finally, for a fermionic pair of phase space variables constructed using a real representa- 
tion of the imaginary unit, with a pair of fermions for each r obeying q r2 = p^ rl , p r -2 = 
but with no relation between p rl and q rl , one simply omits the variables x r2 from d\ip and 
uses the unrestricted form of the fermionic measure for the variables x ri , and the invariance 
argument then proceeds just as before. 
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Appendix B. Canonical Form of {C)av and 
Implications for the Structure of x r and A 

We give here the decomposition of the phase space operators {x r } with respect to 
the canonical form for (C)av given in Eq. (46a) and the relationship this implies between x r 
and x r e ff. We also discuss the canonical form for A which corresponds to that for (C)av ■ 
(1) Real Hilbert space. In real Hilbert space, the anti-self-adjoint operator (C)av is skew 
symmetric, and when the Hilbert space is even dimensional can be brought by a real unitary 
( i.e., orthogonal) transformation to the canonical form {C)av — H ® Cd, with Cd a real 
diagonal matrix and with % 2 the 2x2 skew symmetric matrix 

(When the Hilbert space is odd dimensional the canonical form consists of a block of the 
form i2®Cdi and one further element on the principal diagonal which does not correspond 
to a symplectic structure.) The matrix i 2 spans a two dimensional real Hilbert subspace, 
and a complete set of operators [2] in this subspace can be taken as 1 2 , i 2) W and Wi 2 , with 

from which one sees that i 2 and W anticommute. The operator x r can now be expanded 
over the operator basis provided by Ecjs. (Bl, B2), with coefficients x t a , A = 0,1,2,3 
which are still operators but which commute with i 2 and W, 

x r = x r0 l 2 + x r \i 2 + x r2 W + x r3 Wi 2 . (53a) 

It is convenient to rewrite this expansion in terms of "complex" components denoted by 
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x ra , x rb which both commute with i 2 , according to 

(B3b) 

Writing now = £D, with D nonnegative and with £ a diagonal matrix with elements ±1, 
and writing i e ff = i 2 £ (with a direct product <g> understood), we find 

x r eff = \ [x ra + Wx rb - i2 £(x ra + Wx rb )i 2 £] = \ [x ra + £x ra £ + W(x rb - £x rb £)\ , (54a) 

which reduces in the special case when £ — 1 to 

x r eff = x ra . (BAb) 



(2) Complex Hilbert space. In complex Hilbert space an anti-self-adjoint operator can 
always be written as the c-number i times a self-adjoint operator, and a self-adjoint operator 
can always be brought by a complex unitary transformation to real diagonal form. So we 
have the canonical form (C)av — iCd, writing Cd = £D, with D nonnegative and with £ 
again a diagonal matrix with elements ±1, and writing i e ff — i£ , we have 

x r e ff — -(x r — i£x r i£) = —(x r + £x r £) , (B5a) 

which reduces in the special case when £ — 1 to 

x r e ff = x r . (B5b) 



(3) Quaternionic Hilbert space. In quaternionic Hilbert space the spectral analysis for 
an anti-self-adjoint operator differs in a nontrivial way from that for a self-adjoint operator 
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(see [2] for a detailed discussion and references), and implies that by a quaternion unitary 
transformation, (C)av can be brought to the form ID, with I and D commuting operators 
of the form 

1 =^2\ n } i ( n \ » 



(B6a) 

n\ , 



D=J2\n)D n ( 

n 

and with D n real and nonnegative. We adjoin to this set of operators the additional two 
operators J, K, chosen to commute with D and to form a quaternion algebra with /; one 

possible (but not unique) choice for these operators is 

J = ^2\n)j(n\ , 

" (B6b) 
K = \^\n)k{n\ , 

n 

with i, j, k quaternion scalars. [The J, K of Eq. (B6b) commute with (— 1) F . In a Hilbert 
space with equal numbers of bosonic and fermionic states, so that they can be put into one 
to one correspondence, it is easy to construct an alternative set of operators J, K which 
anticommute with (— l) F .] Let us now expand the operator x r over the basis 1, J, J, K 
with formally real expansion coefficients x t a , A = 0, 1, 2, 3, which commute with the entire 
I , J ,K quaternion algebra (the theory of this is explained in detail in [2]), giving 

X r — X r o + X T \I + X^J + X r %K . (B7a) 

It is convenient to rewrite Eq. (B7a) in terms of formally complex so-called symplectic 
components x ra , x r/ 3 which both commute with /, according to 

{Bib) 

*^ vex — *^ Vo I X f\I j f ^ — Ou ^2 X f^fl • 

Writing now i e ff = I, we get 

x r eff = ^[x ra + Jx r/3 - I(x ra + Jx r/3 )I] = ^[x ra (l + 1) + Jx r/3 (1 - 1)] = x ra . (B8) 
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Although a canonical form for A is not needed for the Ward identity derivation of Sec. 6, 
it is nonetheless instructive to examine the implications that the canonical form for (C)av 
has for A. According to Eq. (45b), there is a function F for which 

(C) AV = F(X, (-If) , (B9a) 

which implies the inverse functional relation 

\=G((C) AV ,(-lf) , 

(B9b) 

=G l ((C) AV , (-If) + (C) AV G 2 ((C) AV , (-If) , 
where G 1 and G 2 are even functions of the argument (C) A y. When we specialize to the form 
(C) A v = ieffD, with D a real positive c-number, as assumed in Eq. (46a), the functions 
Gi i2 reduce to c-number functions of (— \f. Moreover, Gi must be anti-self-adjoint, which 
implies that it must vanish in real and quaternionic Hilbert spaces, giving in these cases the 
canonical form 

\ = i eff D\Gl + G\(-lf] , (510a) 

with G 2 ,x real constants. 

In complex Hilbert space, i is an anti-self-adjoint c-number, so we get the canonical 

form 

A = i[G° + G\(-lf] + i eff D[G° 2 + Gl(-lf] , (5106) 

1 

with G{ 2 re& l constants. To achieve a further simplification in this case, we assume that H 
and H can be expressed in terms of the phase space opertors {x r } using only real number 
coefficients. This implies that (C)av, and hence i e ff, can depend on i only through A, and 
conversely, that A can depend on % only through % e ff\ then the function Gi must vanish, and 
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Eq. (BlOb) reduces to the simpler form given in Eq. (BlOa). More generally, this implies 
that in complex Hilbert space, ensemble averages of monomials constructed from the phase 
space variables can depend on i only through i e ff, a result used at the end of Sec. 6. 



Appendix C. Evaluation of the 5x s Variations of C, H, H, and F 

Varying the definition 

C = ^^x r uj rs x s (Cla) 

r,s 

with respect to x s , we get 

5 Xs C = ^^(5x s u sr x r + x r uj rs 5x s ) , (Clb) 

r 

which by Eq. (10c) becomes 

5 Xs C = u rs (x r 5x s — 5x s e r x r ) . (Clc) 

r 

Varying H with respect to x s and using the definition of operator derivative, we get 

5 x H = Tr—5x s . (C2a) 
ox s 

Applying Eq. (10c) to Eq. (9b) then gives 

S Xa ii = Tr ^ x r u rs 5x s . (C2b) 

r 

We turn our attention next to H = TrH = ReTrif , which is naturally conserved 
under the trace dynamics generated by L = TrL = ReTrL, where L is the same self- 
adjoint operator Lagrangian as appears in the graded total trace Lagrangian L = TrL = 
ReTr(— 1) F L. Since we are assuming that H is conserved under the equations of motion for 
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the graded trace dynamics, the Euler-Lagrange equations for L and L must agree, up to an 

overall sign r\ r = ±1 which can be chosen independently for each degree of freedom, 

5L A 5L _ 

VrPr j 

(C3) 



r . Pr Vr r ■ VrPr ? 

oq r oq r 



5L _ 5L 

Sq r ^ r 5q r 

In fact, for Lagrangians in which fermion time derivative terms all have the structure ip\ ip r , 

no cyclic permutation is involved in varying with respect to ip r since this already stands to 
the right, and so r\ r = 1 for all r and p r = p r , a result that will be assumed henceforth. 
Making a Legendre transformation from L to H, 

H = fr^p r g r -L, (C4a) 

r 

we find 

<5H = Tr J2(<jr$Pr - PMr) ■ (C46) 

r 

Thus the variation of H is given by 

S Xa H = Tr y](q r 6 x ,p r - p r 5 Xs q r ) = Tr(-1) F ^ x r (o rs 8x s , (C5a) 

r r 

with Cj = diag(QB,^B, ■■■,&b) m the notation of Eqs. (10a, b). Forming the sum which is 

needed in the Ward identity derivation, we find 

^ ^ cj us iu) rs =cx ur , 

(C56) 

a=diag(l 2 ,..,l 2 ,-W;..,-W) , 
with 1 2 (for the bosonic variables) and W (for the fermionic variables) as defined in Eq. (B2). 

Carrying out an analog of the discussion of Sec. 3, now using the ungraded trace Hamil- 

tonian H, we find a second conserved operator C given by 

C = ^2x r u rs x s = ^2[q r ,p r ] + ^2[q r ,p r ] . (C6a) 

r,s r,B r,F 
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Defining the auxilliary conserved operator F — — \{C — C), we have 

F = Y,QrPr , (C66) 

which when restricted to the canonical algebra is (up to an additive constant) an anti-self- 
adjoint version of the fermion number operator. To take F into account, we include a term 
TrkF, or equivalently, a term TrAC, in the exponent of the equilibrium distribution; the 
latter makes the algebraic calculations more symmetric and is thus more convenient. We now 
make the essential assumption that the ensemble averages of F and C are functions solely of 
the ensemble average of C and of the grading operator (— 1) F , or equivalently, that the anti- 
self-adjoint operator ensemble parameters k and A are functions solely of A and of (— 1) F , and 
thus commute with A. The validity of this assumption is demonstrated in Appendix G, where 
we study implications of the full Ward identity structure. A straightforward calculation 
shows that 

5 Xs frid = Tr(-l) F [X,J2^rsX r }Sx s , (C7a) 

r 

which has a similar commutator structure to the relation 

S Xs Tr\C = Tr[\,J2 uj rs x r }5x s (C7b) 

r 

deduced from Eq. (Clc). Consequently, the commutativity of the ensemble parameters 
allows us to use the argument of Eqs. (54a-c) to conclude that the variation of the C term 
in the equilibrium distribution does not contribute to the Ward identity. The functional 
relation assumption also implies the continuing validity of the other parts of our analysis 
that depended on the structure of A. Another interesting consequence of commutativity of 
the ensemble parameters is that it implies a vanishing generalized Poisson bracket of the 
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C term in the equilibrium distribution with the C term, with the result that the four first 
integrals appearing in the equilibrium distribution all have vanishing generalized Poisson 
brackets with one another. Inclusion of the C term in the equilibrium distribution of course 
also implies that this term now appears, in a role analogous to that of the C term, in the 
thermodynamic expressions of Eqs. (49a-d). 

Appendix D. Proof of the Lemma of Sec. 6 

We restate and then prove the Lemma used in Sec. 6 to take account of adjointness 
restrictions on the variations. 
Lemma: 

Let Y\ and Y 2 be two self-adjoint bosonic or two anti-self-adjoint bosonic operators 
constructed from the phase space variables. Then in = 5TrY{Y 2 , the self-adjointness re- 
strictions on the variations can be ignored. 



Proof: 

By the cyclic property of Tr, we have 

TrY^ =TrY , 

1 t ( m ) 

Y=-(Y l Y 2 + Y 2 Y l )=Y^ , 

and so it suffices to prove that in = 5TrY with manifestly self-adjoint Y, the self- 
adjointness restrictions on the variations can be ignored. 

We consider first the case of the variation of a bosonic variable x r , for which the self- 
adjointness of x r implies that (Sx r )^ = 5x r . Self-adjointness of Y implies that for each term 
in Y of the form Olx t Or there must be a corresponding term R x r O L , with the grade €l 
oIOl equal to the grade e_R of Or for there to be a nonvanishing graded trace. The summed 
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contribution of the two terms when varied is 

Tve R (0 R L + 0[0 R )5x r , (D2) 

and since the coefficient of 5x r is manifestly self-adjoint we are justified in equating the 
coefficient of 5x r to zero. 

We consider next the variation of a fermionic variable x r , for which there is another 
fermionic variable x s(r ) for which x\ = c r x s ( r ) and x^,-, = c r x r , with c r a c-number of unit 
magnitude with conjugate c r , so that c r c r = c r c r = 1. (The methods for including fermions 
described in Sec. 3 take this form with either c r — — i or c r — ±1.) The corresponding varia- 
tions must thus be related by the self-adjointness restriction 8x\ = c r 5x s ( r y Self-adjointness 
of Y now implies that for each term in Y of the form Olx t Or there must be a corresponding 
term R c r x s (r)0\, with the grade €l of Ol opposite to the grade 6r of Or for there to be a 
nonvanishing graded trace. The summed contribution of the two terms when x r is varied is 
then 

Tr(e R R L 5x r + e L OlO R c r 5x s(r) ) . (D3) 

The self-adjointness restriction on the variations implies that the second term in Eq. (D3) 
is equal to 

Tre L 0[0 R 5xl , (DA) 

which using the fact that Tr of any operator is equal to Tr of the adjoint of the same 
operator, is equal to 

Tr€ L 5x r O R L = Tr{-€ L )0 R L 5x r = Tv€ R R L 8x r , (D5) 

which just doubles the contribution from the first term in Eq. (D3). Hence we get the correct 
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answer by equating the coefficients of 5x r and 5x s ( r ) independently to zero in Eq. (D3). 

Appendix E. Use of the Sources to Generate the Polynomial 
P({x r }) in Eqs. (57d) and (61e) 

We show here that by variation of the sources in Eqs. (57c) and (61c), one obtains 
Eqs. (57d) and (61e), in which J^, 5, [H e ff, ], and [G e ff, ] all act on the Weyl ordered 
polynomial P({x r }) by the Leibnitz product rule. Our argument also applies to Eq. (60b) 
(after dropping the r and r terms), and shows that when V in Eq. (60d) is Weyl ordered, 
the commutator appearing in Eq. (60d) can be evaluated in terms of the canonical algebra 
of Eq. (13) by the Leibnitz product rule. Representing Eq. (57c), Eq. (60b) (with the 
parentheses indicating implicit summation restored and a t replaced by a' t , so that the terms 
inside the expectation read [i e ff(a' t x t e ff),x u ] — Tiuj ut a' t ), and Eq. (61c), after multiplication 
through by the partition function Z, by the generic structure 

= Z(Vx u ) AV , (Ela) 

we wish to show that by varying the sources in Eq. (Ela) we can also derive 

= Z(VP({x r })) AV , (Elb) 

where V acts on the Weyl ordered polynomial P by the Leibnitz product rule 

V(x r x s ) = (Vx r )x s + x r (Vx s ) . (-^Ic) 

We begin by observing that Eq. (Ela) also implies that 

= Z(V(a u x u )) AV , (E2a) 
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and we shall work with the self-adjoint bosonic variables (<J u x u ) henceforth. Multiplying 
Eq. (E2a) by Tr5p u , with 5p u an arbitrary self-adjoint bosonic operator, Eq. (E2a) becomes 

= Z(Tr5p u V(a u x u )) A v , (E2b) 

which we take as the starting point for our discussion. Writing P({x r }) = S[{(a r x r )}], with 
S a totally symmetrized polynomial in its arguments, application of Eq. (Elc) gives 

VS[{(a r x r )}\ = S[{(cr r x r ), r ^ u}; V(a u x u )} , (E2b) 

u 

where the sum over u ranges over the indices of all variables x r which appear as arguments 
of P. Hence to derive Eq. (Elb) it suffices to derive 

= Z{Y,S[{(<J r x r ),r^u};V(<j u x u )}) AV . (E3) 

u 

We begin by varying Eq. (E2b) with respect to the source corresponding to each r^u 
(if a variable x R appears multiple times, we perform multiple independent variations with 
respect to its source pr), and after taking these variations, then summing over all choices 
of u from the among the indices appearing in P. There are two types of source dependence 
which contribute: there are the source terms in the equilibrium distribution of Eq. (51a) 
that we use to form the ensemble averages, and also the explicit source term in the Ward 
identity that was suppressed in writing Eqs. (57c), (60b), (61c), and (E2b). From Eq. (55b) 
we see that the contribution to the summand of this latter term, after left multiplication by 
Tt^Pu^u, is equal to 

s 

with V appropriate to one of the applications discussed in the text. When Eq. (E4a) is 
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varied with respect to the source p v , it contributes 



Z(Tr5p u uj uv a u a,, 



V 



8p v TlY) A V 



(E4b) 



which can be rewritten as 



Z (Tr5 p u 5 p v hV) av^uvO-uO-, 



v 



(EAc) 



Now for each term with the form of Eq. (E4c) in the sum over u, symmetrization implies 
that there is a corresponding term with the roles of u and v interchanged, giving a total 
contribution of 



But this vanishes because for bosonic u, v, the auxilliary quantities a u and a v commute and 
oj uv + oj vu = 0, while for fermionic u, v, we have u uv = u vu and {a u ,a v } = 0. Hence the 
explicit source term in the Ward identity makes no contribution to symmetrized expressions. 
In the remaining terms in the sum over u, the variations 5p u and 5p v each appear in a 
separate graded trace. Thus, after implementing the cancellation of Eq. (E4c), it suffices to 
show that we can derive the generic term in the summand of Eq. (E3) by operations on a 
product of source variation factors, since once the generic term has the correct symmetrized 
polynomial form, all terms in the sum over u are guaranteed to have this form. 

At this point let us take advantage of the fact that the Sp are all arbitrary self-adjoint 
bosonic operators, permitting us to replace them by (— l) F Sp, with the 5p again arbitrary 
self-adjoint bosonic operators, thereby converting all graded traces involving the source 
variations to ungraded traces. We are thus left with the simpler problem of deriving 





= Z{S[{(a r x r ),r ^ u}; V(a u x u )]) AV 



(E5a) 
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given the identity 

= Z(Y[[ReTr{a r x r )5p r }ReTr5p u V{a u x u )) AV . {Ebb) 

To do this, we exploit the arbitrariness of the variations 5p, as follows. Let be a complete 
basis of trace normalized bosonic self-adjoint operators, that is, this basis is characterized 
by the properties that 

TrA A A B = 5 A b , (E6a) 
and that any bosonic self-adjoint operator 0\ can be expanded in the form 

Ox = J2 , 1A = ReTrOiA^ , (EQb) 

A 

which implies the formula 

ReTrdCa = ^ReTrdAAReTrCsA^ . (EQc) 

A 

Now let us take 5pR = {A A , k} and 5ps = A A in Eq. (E5b), with k an arbitrary bosonic 
self-adjoint operator, and sum over A. By Eq. (E6c) this leads to the replacement of the 
product of factors ReTrORSpRReTrOsSps with the single factor 

ReTr{0 R ,K}O s = ReTrK{0 R ,O s } , (E7a) 

which involves the symmetrized (or Jordan) product of the two operators Or^s- Proceeding 
in this fashion, and using the freedom of the k operators just as we use the freedom of the 
operators 5p, we can build up from Eq. (E5b) any trace identity of the form 

= ReTr KZ(S[{(a r x r ),r ^ u}; V(a u x u )\) AV , (E7b) 

with k an arbitrary bosonic self-adjoint operator, which further implies the operator identity 

= Z(S[{(a r x r ),r^u};V(a u x u )]) AV , (E7c) 
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in which S is any self-adjoint polynomial that can be constructed from its arguments by 
repeated applications of the symmetrized product. But repeated application of the identity 

2S(x!, x L ) = S(xi, (x w ), x L )} , (E8) 

w 

with (x w ) indicating that x w is to be omitted from the argument list, shows that any totally 
symmetrized polynomial S can be built up by repeated applications of the symmetrized 
product to its arguments; hence S in Eq. (E7c) can be taken to be S, completing the 
derivation of Eq. (E5a) 

Although we have phrased this derivation entirely in terms of symmetrizing operations, 
it is likely that it can be significantly extended as follows. If we take 5pn = [A a, k], with k 
now anti-self-adjoint, then Eq. (E7a) is replaced by 

ReTr[«, R ]O s = ReTrR[0 R ,O s ] , (E9) 

and thus from the source variations we can in fact build up polynomials which are anti- 
symmetrized in some variables. Moreover, the argument for the vanishing of the explicit 
source term contribution requires not total symmetrization, but only symmetrization in all 
aguments x u , x v for which the symplectic structure oj uv is nonzero. Thus, if we define a 
partially Weyl ordered polynomial to be a polynomial which is symmetrized with respect to 
all arguments x u , x v for which u uv ^ 0, then it appears likely that with careful attention 
to the combinatorics, one should be able to show that Eq. (Ela) implies the extension of 
Eq. (Elb) in which P is any partially Weyl ordered polynomial in its arguments. The ex- 
ample of Eq. (17e) also indicates that there will be extensions of Eq. (Elb) to certain cases 
in which P is a non-Weyl ordered polynomial of low degree. We expect these cases to play 
an important role in the study of operator gauge invariant theories, which we will take up 
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in detail elsewhere. 



Appendix F. Canonical and Symmetry Transformations 



We give here further details of the structure of canonical transformations in generalized 
quantum dynamics, with special emphasis on their role as symmetry transformations. 

A particularly interesting class of canonical transformations are what we shall term 
linear symmetry transformations, defined as transformations of the type Eq. (27) generated 
by total trace functionals of the special form 



These transformations linearly transform the canonical coordinates {q r } among themselves, 
with a corresponding transformation on the canonical momenta, but do not mix coordinates 
with momenta. When the indices r, s in Eq. (Fla) are both bosonic or both fermionic, the 
coefficients h rs are taken to be ordinary c-numbers, while when one index is fermionic and 
one is bosonic, the coefficients h rs are taken to be Grassmann c-numbers. Thus, the linear 
symmetry transformations include grade-changing transformations which mix the bosonic 
and fermionic coordinates. As a consequence of the grading structure of h, we have 



which together with the cyclic property of Tr implies that Gh is Weyl ordered. Thus, a linear 
symmetry transformation is also a Weyl ordered intrinsic canonical transformation. Under 
the generalized Poisson bracket operation, two linear symmetry transformations compose as 
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G h =TrG h 




(Fla) 



r,s 




(Fib) 



[1,2] 



{Gg, G h } - G[g jh ] 



(F2a) 



with [g, h] the commutator 



[g, h] rs = ^2(g r th ts - Ktgts) 



(F2b) 



and hence linear symmetry transformations form a Lie commutator algebra under the gen- 
eralized Poisson bracket. Clearly a linear rearrangement of the canonical coordinates among 
themselves, with c-number or Grassmann c-number coefficients, together with a correspond- 
ing linear transformation among the momenta, transforms a Weyl ordered polynomial into 
another such polynomial. Therefore the set of Weyl ordered total trace functionals, and 
thus of Weyl ordered intrinsic canonical transformations, is closed under the action of linear 
symmetry transformations. [We note in passing that in Refs. [1, 2] we also introduced linear 
symmetry transformations in which h is an arbitrary quaternionic (hence non-commutative) 
coefficient matrix; although the Lie property of Eqs. (F2a, b) holds for this generalization, 
most of the other properties of canonical transformations derived in Sec. 4 and this Ap- 
pendix do not. For example, symmetry transformations based on quaternionic representions 
of compact Lie groups do not leave C invariant.] One can also define a generalization of 
linear symmetry transformations with generators which can mix coordinates and momenta 
according to 



formations are also Weyl ordered. These transformations also form a Lie algebra under the 
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G h =TrG h 




(F2c) 



r,s 



with the grading structure of h rs analogous to that for h rs , which implies that these trans- 



generalized Poissson bracket, with the structure 



{ G 9' G h} ~ G k 



hu = 2,{gtr + t r g rt )uj rs (h 



(F2d) 



r,s 



certain Bogoliubov transformations are of this more general type. 

We shall now derive a second relation which is similar in structure to Eq. (29) of the 
text, and which describes the action of a linear symmetry generator = Tr Gh on a Weyl 
ordered intrinsic canonical generator G = TrG, when the phase space operators {x r } are 
specialized to the canonical algebra, 



Since by the Weyl ordering hypothesis G is symmetrized, we can represent it as a sum of 
monomial terms produced by generating functions with the form g n of Eq. (11a), and it then 
suffices to prove the identity for only one such term. Writing g in the form 



{G,,G} 



Tri[G h ,G] 



(F3) 




(F4) 



r 



use of the generalized Poisson bracket in the form given in Eq. (6b) gives 



{G h ,Trg n }=nTrg h g 



n-l 



(F5a) 



with g h defined by 




(F5b) 



r,s 



One can now check that over the canonical algebra one has 




(F5c) 



r,s 
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But the commutator [— iGh, g n ] reduces to the totally symmetrized product of gh with <? n_1 , 
which is equal to the right hand side of Eq. (F5a) under the trace, completing the proof 
of Eq. (F3). Thus, there is an isomorphism between (a) the action of a linear symmetry 
transformation on an arbitrary Weyl ordered intrinsic canonical transformation, under the 
generalized Poisson bracket operation of generalized quantum dynamics, and (b) the corre- 
sponding behavior of the canonical algebra specializations of these transformations, under 
the usual commutator operation. This isomorphism extends to the more general Bogoliubov 

type transformation of Eq. (F2c), where we find 

{G h ,G} = -Tr i [G k ,G]=nTrg k g n - 1 , 

~ ~ (F6) 

9h = 2_^(uj rt a t h rs x s + x r h rs u st a t ) . 

r,s,t 

The isomorphism does not extend, however, to the action of generic Weyl ordered 

intrinsic canonical transformations on one another. To see this, let us consider the case 

of two such canonical transformations with generators Gi and G 2 which are generating 

functions for Weyl ordered monomials, 

d = Trd, G 2 = TrG 2 , 



Gi — <?™\ G 2 — 2 > 

<?1,2 = ^]cTl,2r x r ■ 



(F7a) 



Then from Eq. (9a) we find 



{Gi,G 2 } =Tr^n 1 ^ 1 l a lr LU rs n 2 g% 2 1 a 2s 



r,s 

— 1 r ,ri2 — 1 



(F7b) 



=CTrg?- l g? 
C =nin 2 ^^cr lr uj rs a2s , 

r,s 

which is clearly not Weyl ordered when ri\ and n 2 are both greater than 2. Thus the 
generalized Poisson bracket of the generators for two Weyl ordered intrinsic canonical trans- 
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formations is in general not a Weyl ordered canonical transformation. On the other hand, if 
we specialize to the canonical algebra and then evaluate the commutator of G\ and G 2 , we 
easily find using Eq. (13) that 

[Gi, G 2 ] — iCS(g'^ 1 ~ 1 , #2 2_1 ) , (F7c) 

with S the polynomial formed from completely symmetrizing n\ — l factors of g\ with respect 
to n 2 — 1 factors of g 2 , which is Weyl ordered. In other words, over the canonical algebra, the 
commutator of two Weyl ordered generators is Weyl ordered, but this does not correspond 
to the composition properties of Weyl ordered canonical generators under the generalized 
Poisson bracket operation. 

In future work we plan to give a more detailed analysis of the Poincare transformations 
in generalized quantum dynamics than was given in [1, 2], including a study of their rela- 
tionship to the classification of canonical transformations given in Sec. 4 and here. 



Appendix G. The Full Ward Identity Structure 

We describe here the full Ward identity structure resulting when the presence of the ad- 
ditional conserved anti-self-adjoint operator F of Eq. (C6b), or equivalently C of Eq. (C6a), 
is taken into account. Including the latter in the equilibrium distribution, Eq. (51a) becomes 

p =Z~ X exp[-Tr ^p r [o r x r )\ exp(-frA(7 - TrAC - f H - rH) , 
Z= dpexp[-Tr^p r (a r x r )]exp(-fr\C -Tr\C -tH-tH) ; 

J r 

ensemble averages denoted by the subscript U AV" will be understood henceforth to be taken 
in the distribution of Eq. (Gl). Including the C term in the Ward identity of Eq. (55b), and 
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for the moment retaining all terms coming from the variation of the equilibrium distribution, 

we get 

=(([\,x u ] + (-1) F [A, ^a ur x r ] + r(-l) F ^ a ur x r + tx u + ^u us a s p s )TiV 
1 5V 

s 

giving the full version of the Ward identity derived in Sec. 6. Let us now take V = H, 
and make the approximation (expected to be valid in the large N limit) of replacing H by 
its ensemble epectation (H) AV . The coefficient of this term in Eq. (G2) then becomes the 
proportional to the variation ^ s u us 5Z/5x s of the partition function of Eq. (Gl), which is 
zero, and so Eq. (G2) reduces, after use of the equation of motion of Eq. (9b) and division 
by h, to 

^ = {{]^{h~ l ieff,H},x u }- x u )av ■ (G3) 

This is the same relation as was obtained in Eq. (56b) of the text, but we have now shown 
that its derivation does not depend on the assumption that the ensemble parameters A and 
A are functionally related. 

Let us now derive a second Ward identity by using properties of the operator C. Since 
this is an anti-self-adjoint operator, its ensemble expectation can be written in a polar form 
analogous to Eq. (46a) for (C)av, 

(C)aV = ieffD , ieff = -iff , ilff = ~1 , 

(G4) 

[i e ff, D] — , D self — adjoint and nonnegative. 

The difference is that we do not now assume D to be diagonal, but we shall assume it to be 
nonsingular, so that the inverse (-D) -1 exists. Let us now derive a second Ward identity, by 
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considering the expression 

= J dfi5 Xs [pTr{d, (D)- i: i eff }vj , (G5a) 

with p the full equilibrium distribution of Eq. (Gl). Proceeding as in Eqs. (53-55) of the 
text, with the one exception that we now multiply through by \ J2 S u us , we end up with the 
following analog of Eq. (G2), 

=(([A, X u \ + (-1) F [A, ^ a urX r ] + t(-1) F ^ a ur%r + TX U + (-if ^ ^us^sPs)^ 

r 

+[k(D)-% f ,V},x u ]-Y,i 



,Wusj — )AV , 



(G56) 

where the hatted operator derivative SY/Sx s is defined by 

5V = frJ2^^ s . (G5c) 

8 5X » 

Let us now take V — H in Eq. (G5b), and replace the conserved extensive quantity 
H by its ensemble average (H)av- The coefficient multiplying this quantity in Eq. (G5b) 
now becomes proportional to the partition function variation ^ s u us 5Z/Sx s , which again 
vanishes, and so Eq. (G5c) reduces in this special case to 

1 

0=(L{(D)- 1 l ff ,H},x u ]-J2^^)^ ■ ( G6fl ) 
z „ ox s 

s a 

But comparing the definition of Eq. (G5c) with Eq. (C4b), we see that the fact that the 
same equations of motion follow from the ungraded and graded trace variational principles 
can be expressed succinctly by the identity 



. ^ 5H ~ 5H 

/ —~ ox s <w 

s s s 
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and so Eq. (G6a) simplifies further to 

= ([\{(D)-%f, H}, x u \ - x u ) AV . (G7a) 

Comparing this expression for the ensemble average of x u (still in the presence of sources!) 
with the similar expression obtained from the original Ward identity in Eq. (G3), we conclude 
that we must have the relation 

n-\ Jf = {b)-% J , {Gib) 

which since D is nonnegative implies the further relations 

D = fi, i ef f = i ef f , 

(G7c) 

(C) AV = (C) AV , (F)av = , 
justifying the functional relation assumption made in Appendix C. Substituting Eq. (G7b) 

into the second Ward identity of Eq. (G5b), and multiplying through by H, we get 

=(([A, X u ] + (-1) F [A, ^ a urX r ] + t(-1) F ^ a ^r + TX U + (~l) F ^ ^us^sps^V 

r r s 

+ [li i eff,V},X u ] -h^2LJ us ^-) AV , 

s s 

(G8) 

giving the form of the second Ward identity analogous to Eq. (G2). 

We must now perform an important consistency check. If instead of taking V = H in 
the Ward identities of Eqs. (G2) and (G8), we take instead V = (—1) F H, the effect is to 
simply interchange the roles of H and H in the argument showing that the variation of the 
equilibrium distribution does not contribute, and so this argument remains valid. Hence we 
get two new relations, which we must check are not in contradiction with Eq. (G3). The two 
new relations are seen to be identical when one uses the identity 

<JX S te, 



s s 
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which like Eq. (G6b) is a consequence of the fact that the same equations of motion follow 
from the graded and ungraded trace variational principles. When the index u is bosonic, 
the new relation is easily seen to be identical in form to Eq. (G3), and so is automatically 
consistent. When the index u is fermionic, we introduce the definition H e ff as in Eq. (57a), 
and rewrite Eq. (G3) in the form 

{[H e ff,x u ])Av = h(^2u us — ) AV ; (G9b) 



s s 



in this notation, the new relation takes the form 



{{H ef f,x u }) A v = hi^Lbus — )av , (C9c) 



Sx s 

s 



which involves a commutator rather than an anticommutator on the left hand side. Treating 

separately the fermionic cases in which u — 2r — 1, x u = q r and u = 2r, x u = p r , and taking 

linear combinations of Eqs. (G9b) and (G9c) which eliminate the right hand side, we find 

that the new relation of Eq. (G9c) implies that 

{H e ffq r )Av =0 , 

(G10) 

(p r H ef f) A v =0 . 

These relations are compatible with what one expects for vacuum expectations in quantum 
field theory, when one defines the fermionic vacuum so that (0\p r e ff = Qr e//|0) = 0; hence 
the new relation of Eq. (G9c) is consistent with the isomorphism conjectured in the text. 
We also remark that Eq. (G7c) for the expectation of F, which can be rewritten as 

C}2qrPr)AV = , (Gila) 

r,F 

and Eq. (G10) do not contradict the fermionic canonical algebra of Eq. (60c), because in 
general 

q r p r ^ q r e ffp r e ff ; (Glib) 
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referring to the analysis of Appendix B we see that equality in Eq. (Glib) can hold only for 
the special case [see Eq. (B5b)] of complex Hilbert space with £ — 1, which is thus ruled 
out. 

Having established the equality of the ensemble averages of C and C, and therefore a 
functional relationship between the corresponding ensemble parameters A and A, and (— 1) F , 
we can apply the argument of Eq. (54c) to conclude that the commutator terms involving 
A and A drop out of the Ward identities. The Ward identities of Eqs. (G2) ad (G8) then 
simplify to 

={{i-(-l) F ^a ur x r + rx u + ^2u us a s p s )hV 

i r ^ sv ( G12fl ) 



and 



=((r(-l) F ^ a urXr + TX U + (~1) F ^ U us (T s p s )hV 

1 5V {GUb) 

+ [^{ i eff,V} ) X u ] -h^2uJ us j—) AV , 

s " X s 

with Eq. (G12a) the basis of further applications as discussed in the text. For any V obeying 
the conditions 

x ^ SV _ ^ A 5V 

(GUa) 

s s a 

an analysis paralleling the consistency check on the Hamiltonian given above shows that 
the Ward identities obtained from Eqs. (G12a, b) when the f and r terms are neglected, 
and those similarly obtained when V is replaced by (— 1) F V, are consistent as long as the 
conditions 

(V eff q r ) AV =0 , 

(G136) 

(PrVeff)AV =0 , 
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with V e ff = \{i e ff,V}, are obeyed for fermionic r. These conditions are compatible with 
those of Eq. (G10), and have the same interpretation in terms of fermionic vacuum structure. 
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